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If the Electroweak Symmetry Breaking Sector turns out to be strongly interacting, the actively 
investigated effective theory for longitudinal gauge bosons plus Higgs can be efficiently extended 
to cover the regime of saturation of unitarity (where the perturbative expansion breaks down). 

This is achieved by dispersion relations, whose subtraction constants and left cut contribution can 
be approximately obtained in different ways giving rise to different unitarization procedures. We 
illustrate the ideas with the Inverse Amplitude Method, one version of the N/D method and another 
improved version of the K-matrix. In the three cases we get partial waves which are unitary, 
analytical with the proper left and right cuts and in some cases poles in the second Riemann sheet 
that can be understood as dynamically generated resonances. In addition they reproduce at Next 
to Leading Order (NLO) the perturbative expansion for the five partial waves not vanishing (up 
to J=2) and they are renormalization scale (/x) independent. Also the unitarization formalisms are 
extended to the coupled channel case. Then we apply the results to the elastic scattering amplitude 
for the longitudinal components of the gauge bosons V — W, Z oX high energy. We also compute 
hh hh and the inelastic process VV hh which are coupled to the elastic VV channel for 
custodial isospin 7 = 0. We numerically compare the three methods for various values of the low- 
energy couplings and explain the reasons for the differences found in the I — J — 1 partial wave. 

Then we study the resonances appearing in the different elastic and coupled channels in terms of 
the effective Lagrangian parameters. 

PACS numbers: 12.60.-i, 12.60.Fr, 12.39.Fe 


I. INTRODUCTION 

The most outstanding discovery in particle physics during the last years is probably the finding that the LHC 
collaborations ATLAS [T] and CMS [2) published in 2012 announcing a new boson with scalar quantum numbers and 
couplings compatible with those of a Standard Model Higgs at about 125 GeV. However the first LHC run finished 
without any other finding [3] up to an energy of 600-700 GeV (and higher yet for additional vector bosons). This 
lightness respect to any new physics could alternatively suggest that the Higgs is indeed an additional Goldstone boson 
(together with those giving rise to the W and Z masses) related with some global spontaneous symmetry breaking 
triggering the SU (2) p x U{l)v U{l)em gauge symmetry breaking [4]. 

If that were the case, some effective description of the Symmetry Breaking Sector (SBS) of the Standard Model 
(SM) would be appropriate (see for example [5IIT5]). The presence of that energy gap is also suggestive of a non-linear 
realization (the most general approach to the effective theory). The old Electroweak Chiral Lagrangian (ECL) |I6| 
technique, based on standard Chiral Perturbation Theory (ChPT) for QCD [17], can be extended to include the new 
found Higgs-like light particle as a scalar singlet. 

In a recent work m we have shown that, for essentially any parameter choice except that of the Standard Model 
and perhaps other very carefully tuned sets, the interactions will generically become strong at sufficiently high energy, 
and have argued that a second, very broad scalar pole is expected. In a more recent work m we performed the 
one-loop computation of the two-body scattering amplitudes among the uj Goldstone bosons and the h scalar by using 
a generic effective Lagrangian, in the kinematic regime M| s < 47ru ^ 3 TeV. 

In spite of the success of one-loop computations (for example in ordinary ChPT) it is clear that it can only be useful 
at very low energies. Moreover, the case in point of the ECL deals with the would-be Goldstone bosons uj that are 
eventually related to the longitudinal components of the gauge bosons only through the Equivalence Theorem (ET) 
123 , which is valid only in the kinematic regime M| s, that corresponds to a high energy limit. 

Thus even at low energies one could not expect the truncated series to apply in a context of strong interactions. This 
situation is not improved in any significant way by going further in the chiral expansion by computing two or more 
loops. Going to higher orders one has to deal with a very fast increasing number of chiral couplings and extremely 
complicated computations. 

On the other hand one can try a different strategy to extend the low energy regime by using dispersion relations 
(DR) compatible with analyticity and unitarity. This program has proved to be extremely useful in the original 
ChPT as applied to low energy hadron dynamics and it is quite possible that it could also be useful for the SBS of 
the strongly interacting sector of the SM if properly applied. This program is bringing much attention, because some 
predictions can be checked at the LHC Run-II [5], as we will show here for the one-loop computation. However, the 
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use of unitarization methods for extending the applicability of ChPT has been criticized because they are considered 
arbitrary in some sense, and the results they provide depend on the particular method considered. 

In this work we will try to show that the one-loop results, when properly complemented with DR, can provide an 
analytical and unitary description of higher energy dynamics which is essentially unique qualitatively; at least so up 
to the first resonances which can also be described as poles in second Riemann sheet due to the proper analytical 
behaviour of the amplitudes. The rationale for this is that any physically sensible amplitude must fulfill the appropriate 
DRs which are typically integral equations. In principle those equations have many solutions. However one can impose 
some particular dynamics by performing subtractions on the DR and fixing the values of the subtraction constants. 
In our case it is clear that these constants must be obtained from the effective Lagrangian and so we will introduce 
the different dynamics compatible with the low-energy expansion in the DR relation with an appropriate number of 
subtractions. In this way one expects to reduce enormously the space of possible solutions of the integral equations, 
at least up to the first poles. If this is really the case, different unitarization methods will provide qualitatively similar 
results and the differences could be attributed to the different approximations used for solving the DR equations. 

We study in detail the Inverse Amplitude Method (lAM, formerly called the Fade method) [HI [22], the N/D [23] 
method as applied to the SBS of the SM and also the so called improved K-matrix method (see [24],[25| and [26| 
for exposition and some uses of the non-improved method). The main novelty in the case of the electroweak chiral 
descriptions compared with standard ChPT applied to hadrons is that now the GB are really massless. We then pay 
attention, not only to the usual ultraviolet (UV) divergences appearing in the DR integrals but also to the infrared 
(IR) ones. 

Thus we rederive the lAM method for massless particles using a twice-subtracted DR instead of the original 
derivation that used three subtractions [22] • For our construction of the N/D method we introduce a renormalization- 
scale invariant splitting of the NLO amplitude into a left and a right parts each bearing the corresponding left and 
right cuts. Then we write a thrice-subtracted DR for the denominator function and solve the corresponding integral 
equation by iteration (in fact one is sufficient to get a sensible result). Finally we consider also an improved version of 
the K-matrix method [ 111127 ] that produces partial waves having a proper RC that allows for analytic continuation 
to the second Riemann sheet in the search for poles (resonances), thus fixing the typical absence of these poles of this 
unitarization method. 

We follow the natural order of presentation, with the Effective Lagrangian briefly recounted in subsec. [TTA| followed 
by a short discussion on the elastic and inelastic scattering amplitudes in subsec. m A part of this work can be 
found in our earlier article[T9]. We have however calculated the fifth, non-vanishing NLO amplitude with angular 
momentum and custodial isospin 2, a new result not commonly quoted in analogous hadron systems. Thus, we have 
now exhausted the massless low-energy Lagrangian to NLO, by computing all nonvanishing channels and interchannel 
couplings. The calculated amplitude coefficients, their behavior under scale changes, and our conventions for the 
partial wave amplitudes, are all given in the appendix to make this subsection more readable. 

We dedicate section [Till to the Inverse Amplitude Method, both for single and coupled channels, and provide a new 
derivation based on twice-subtracted dispersion relations especially useful for massless particles. Section [TV| in turn is 
dedicated to generically describing the N/D method, but also to construct a new approximate solution thereof that 
has the same desirable physical properties of the lAM (and unsurprisingly, both coincide where both are applicable). 
More so, we explore the Improved-K matrix method in subsection |V A| and compare all three methods extensively. 

The computer evaluations of all three methods are exposed in sections |V^ and m and we conclude that the 
unitarization methods are qualitatively robust and a reliable guide in the search for strongly-interacting new physics, 
with little model dependence. 

After a terse summary in section [VlTTl we dedicate four appendices to technical details of the perturbative amplitude 
calculations, their partial waves, the coupled-channel lAM, the (one-iteration) solution of the N/D method, and the 
numerical extraction of poles in the complex plane (resonances if in the second Riemann sheet or tachyons if in the 
first Riemann sheet). 


II. ELECTROWEAK CHIRAL LAGRANGIAN AND SCATTERING AMPLITUDES WITH A HIGGS 


We have already presented the Lagrangian density and perturbative LO amplitudes in m for VV and hh scattering. 
Here we quickly remind of the basic equations with reduced discussion and settle to a more standard notation than 
we previously used. Also in m we obtained the one-loop scattering amplitudes between the (massless) uj and 
h. This section is divided in two subsections, one |H A| dedicated to exposing the effective Lagrangian and another 
subsection H B| dealing with the scattering amplitudes. Some further detail is relegated to appendix [A| 
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A. Effective Lagrangian 


There are several equivalent forms of the universal Electroweak Chiral Lagrangian employing only the experimentally 
known particles. At leading order we adopt the gauged SU{2)l x SU{2)r/SU{2)c = SU{2) Non-linear Sigma 
Model (NLSM) coupled to a scalar field h as 

Co = + ^d^hd>^h - V{h) ( 1 ) 

with U = + iuj/v] uj = uJaT^ parameterizing the would-be Goldstone boson (WBGB) field. Since we will 

neglect the coupling to transverse gauge bosons, ^ in this article’s computations. 

The constant v is well known from Fermi’s weak constant, := 1/{\/2Gf) = (246 GeV)^. The scalar field interacts 
through T, an arbitrary analytical function; in the effective-theory approach we need only the first terms of its Taylor 
expansion 

T{h) = l + 2a-+b(-'] +... (2) 

V V'^/ 

widely used in the literature^. 

Reference [28] provides some recent experimental bounds on a and b that we employed in [18]. Finally V is an 
arbitrary analytical potential for the scalar field, that is of no further reference in this work. 


V{h) = Vnh^ = Vo + 

n=0 


1 




M 2 


h 1.3 


Mi 


h 


2v 


Sv‘^ 


At NLO in the chiral expansion we need to add the four-derivative terms 
£4 = a^itrV^V.f + 

+ ^{d^hd^hf + ^{d^hd>^h)tr{D^U)^D''U + A,{d^hd''h)tr{D^U)^D^U - 


(3) 


(4) 


where = D^UW. We have explicitly written only the five terms strictly needed to renormalize the one-loop elastic 
WBGB scattering amplitudes (for s ^ and the coupled-channel processes uouj hh and hh hh. These terms 

produce additional contributions to the amplitudes which are of order s^. 

The a 4 and chiral parameters multiply the operators Odi and Od 2 in the classification of m Those, as well as 
the additional ones d and e, encode the dependence on the possible underlying dynamics triggering the spontaneous 
symmetry breaking of electroweak interactions. They all vanish in the MSM. The operators with coefficient d and e 
correspond to Oi and O 2 as classified by Azatov et al. [7], or to P 19 and P 20 in [9], and are NLO equivalent to Odt, 

Od8 in mi 

The operator associated with g is denoted as Ph in |H] and Odii in mi In line with recent literature [ansi , we 
are not much interested in the dynamics of the pure Higgs sector, since the process hh hh will hardly be measured 
in the foreseeable future. But since the NLO renormalization of our effective Lagrangian requires this one operator, 
we will assess its numeric effect on the ujuj channel in figure below, where we see that for it to be sizeable the values 

of g have to be quite unnaturally large. _ 

We have given the renormalization of these operators in na, and we rewrite it in the new notation in appendix | A 2 
An off-shell analysis that covers a larger number of operators has also recently appeared m- 

The Lagrangian in Eq. 0 and Eq. 0 contains the more general low-energy physics of the ESBS for any conceivable 

dynamics having at least an approximate SU ( 2 ) custodial isospin symmetry in the limit g = g' = 0. 

The easiest example is the Minimal SM (MSM) [31], which corresponds to the parameter selection a = b = 1 and 
a 4 = as = ^ = d = e = 0. The Higgs field H is just the scalar field h, so that Mjj = = 2Au^, and the scalar 

self-couplings are ds = Au, ^4 = A/4 (both proportional to M|) and = 0 for i > 4. 


^ One can alternatively employ T{h) = +---7 where / is an arbitrary, new-physics energy scale, as we have done in 

recent work. This is perhaps more natural if the Higgs happens to be the Goldstone boson of a higher symmetry broken at the scale 
/, but in this article we adopt the more widely used convention of employing v, the SM symmetry-breaking scale. Obviously a = olv/ f 
and b = I• 

^ The two operators multiplying d and e are of dimension 6 in what would concern transverse gauge-boson inelastic scattering WtWt —>■ hh, 
but are of dimension 8 for the longitudinal ones, as seen upon expanding U as we will show shortly in Eq. 
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Another interesting class of models are the dilaton models [31] where h would represent the dilaton field and we 
have o? = h = v‘^/p with / being the scale of the scale symmetry breaking. The potential and NLO parameters 
depend on the particular dilaton model but in any case the di are also of order for any i. 

Third, the popular SO(b)/SO{A) Minimally Composite Higgs Model [33] (MCHM) also provides an example where 
o? = 1 — v‘^ /P and 6 = 1 — 2v‘^/p (where / is in this case the scale of the 5'0(5)/5'0(4) symmetry breaking) while 
the scalar-boson self-couplings contingent on model details, are of order too. 

Finally it is also possible to reproduce the old Higgsless Electroweak Chiral Lagrangian (EWChL) in [16] by the 
simple parameter choice a = b = g = d = e = 0. 

To address the high-energy (i. e., for Ps ^ 100 GeV) elastic scattering of the longitudinal components of the 
electroweak bosons, we can apply the ET m, which relates the WBGB amplitudes with the corresponding amplitudes 
involving longitudinal components on the electroweak bosons at high energies. Eor example one has: 

T{WlWl^WlWi)=T{u'^u’’+ , (5) 

Thus the ET allows to carry the computations out with the simpler WBGB dynamics. This theorem applies for any 
renormalizable gauge, but for the Landau gauge (where the WBGB are formally massless) it is especially useful and 
transparent. Since the transverse degrees of freedom are weakly coupled to the longitudinal sector, to explore just 
the latter we will set g = g' = 0. 

The remaining active degrees of freedom are the massless (Landau-gauge) WBGB, and the Higgs-like scalar h that 
will be considered massless in the following as we are interested in the high energy region. According to ATLAS and 
CMS, Mh — 125 GeV. Then ^ Mw ^ Mz ^ 100 GeV and consistency requires to consider the massless h limit, 
i.e., Mh — 0 if one is only interested in the energy region where the ET can be applied. Consequently we concentrate 
on WBGB scattering for M^, M^, Mf 0 <C s < A^ where A is some ultraviolet (UV) cutoff of about 3 TeV, setting 
the limits of applicability of the effective theory. 

According to the results of LHG Run-I I3E1], no new physics has been discovered up to an energy of about 600- 
700 GeV. However, the center-of-mass energy of the LHC is going to be increased from 7-8 TeV to 13-14 TeV (and 
the luminosity will be much higher) at Run-H. Thus, the applicability limit s < (3TeV)^ of the theory is 

within the new range of energy. Actually, LHC Run-H is going to be a great opportunity to check strongly interacting 
EWSBS theories controlled by unitarity [5]. 

We will also assume that the di self-potential parameters of the h-scalar are of order so that we can neglect 
them altogether, as is natural in the three particular models just mentioned, MSM, Dilaton, and CHM. 

Under these kinematics, Eqs. 0 through 0 yield a Lagrangian density 

+ + ^(d.hd^hf 

2d 2 p 

+ ^d^hd^hd^uj^d^uj^ + ^dpd^hd^uj^d^uj^. ( 6 ) 


B. The WBGB scattering amplitude in EWChPT at the one-loop level 

Concentrating first on elastic scattering, the custodial symmetry of the SBS of the SM in the limit g = g' = 0 
allows to write the WBGB amplitude uJaUJi) ujc^d as 

'^abcd ~ L P^ah^cd 4“ ^(C <^5 P^aPbd 4“ A('U, t, Sp^d^bc • P) 

Because of crossing symmetry for four identical particles, only one function of the Mandelstam variables A is needed. 
In terms of the charge states uj^ = {P ^ iP)lp2 and z = P the amplitudes can be written as: 

A{uj~^uj~ ^ zz) = A{sA^u) (8) 

A{uj~^uj~ ^ = A(s, t, rt) + A(t, s, 1^) 

A(zzzz) = A(s,t,p + A(t, s,u)-j-A(u,t, p 

(the remaining charge combinations can be obtained from these by crossing symmetry). The ApA^u) amplitude can 
be expanded in a similar way to ordinary ChPT. Quoting the NLO tree-level and one-loop subamplitudes yields 

^ = A( 0 ) + ^( 1 ) . . . = ^( 0 ) + + ^( 1 )^ , , , 


(9) 
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The next two-body processes to consider are the channel coupling uja^b hh between two uj WBGB and a scalar 
boson pair and hh that are needed to obtain one-loop unitarity in ujuj scattering. Obviously both processes 

have the same amplitude because of time reversal invariance. With h being an isospin singlet, the amplitude takes 
the form 


Mab{s,t,u) = M{s,t,u)Sab 

We also consider the amplitude for elastic scattering hh hh^ 

r{s, t, u) = 


■-tree + ^1 


( 1 ) 

loop ' ■ ' 


( 10 ) 


( 11 ) 


All these amplitudes are explicitly given in appendix ] A 1[ 

The unitarity of these three scattering amplitudes is best exposed in terms of the isospin- and spin-projected partial 
waves; this requires projecting over custodial-isospin and angular momentum. For elastic WBGB scattering there are 
three custodial-isospin Aj amplitudes (/ = 0,1, 2), analogous to those in pion-pion scattering in hadron physics, 


Aq{sA^u) = 3A(s, t, 1 ^) + A(t,s, 1 ^)-h A(i 4, t,s) 

Ai{sA^u) = A{t^s^u) — A{uA^s) 

A2{sA^u) = A(t,5,rt) + A(rt,t,s) . 

The projection over definite orbital angular momentum (the WBGBs carry zero spin) is then 


( 12 ) 


Auis) = —— 


■e/> 


cosO) Pj{cos 0) Ai{s A • 


These partial waves also accept a chiral expansion 

Auis) = Afj{s)+A<-^]{s) + ..., 


(13) 


(14) 


that take the general form 


Af]{s) = Ks 

4V(s) = (B{p) + D\og^ + ElogA 

V fi fi 


(15) 


The constants K, D and E and the function B{jii) depend on the different channels IJ = 00,11, 20, 02, 22, as shown 
below in appendix ] A 1[ and we will use the same notation for the inelastic and pure-h scattering reactions. 

As Aij{s) must be scale independent we have 


B{ii) — B{fio) P {D P E) log ; 

ho 


(16) 


This B{jii) funct ion d epends on the NLO chiral constants (with certain proportionality coefficients p 4 and ps that can 
be read off Eq. (A17) and following) 


B{n) = Bo + p4a4{ij) + p^aziij) , 


(17) 


where Bq also depends on a and b and from now on we omit the superindices r on the renormalised coupling constants 
for simplicity. 

Since the “Higgs” boson is assigned zero custodial isospin, ujuj hh and hh hh occur only in the isospin zero 
channel / = 0. 

The normalisation of the \ujuj)i=o state introduces a factor 1/v^ and the sum over the three contributing charge 
combinations (H—, —h, 00) a factor 3, so that for the inelastic amplitude we have Mq{ujuj hh) = \/3M{s, t, u). For 
the scalar-scalar interaction there is no such factor and To{hh hh) = T{sA^u). Omitting the isospin subindices 
(which take only the value 0) and proceeding to the angular momentum projections, we find the chiral expansions 
equivalent to the ujuj elastic one in Eq. (15). They read 


Mj{s) = K's+{B'{pi)+D'\og^+E'\og^]s‘^ + ... 


Tj{s) = K" 




B"{p) + D"\og^+E"\og^ 


(18) 
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(with J subindex omitted in the constants). The functions and are in all analogous to B{jii) in Eq. 

renormalization is carried out by d, e (for B') and g (for B") involving the h boson. 

The partial-wave amplitudes Ajj{s), Mj{s) and Tj{s) are all analytical functions of complex Mandelstam-s, having 
the proper left and right (or unitarity) cuts, shortened to LC and RC respectively. The physical values of their 
argument are s = + ie (i.e. on the upper lip of the RC), where Eqm is the total energy in the center of mass 

frame. For these physical s values, exact unitarity requires a set of non-trivial relations between the different partial 
waves that we now spell out. For the problem of ujuj scattering considered here the reaction matrix is block-diagonal: 


F{s) 


/Eoo 

0 

0 

0 


0 

Eo2 

0 

0 

0 

0 

0 

En 

0 

0 

0 

0 

0 

^20 

0 

\ 0 

0 

0 

0 

J 


where Fjj{s) are the partial-waves matrices. For example for 1 = 0 we have; 


Poo{s) 


f ^ 00 (s) Mo{s)\ 
\Mo{s) Ms)) 


and 


Fo2{s) 


(AMs) M 2 {s)\ 
[M2{s) Ms))- 


For / 7 ^ 0 there is no mixing with the hh channel and the Fij{s) matrices have just one single element: 


Fij{s) = Ajj{s) 

Now unitarity requires that on the right cut: 

ImF(s) = F{s)F\s). 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 

(23) 


This equation produces a set of relations concerning the different partial waves. For 1 = 0 and either J = 0 or J = 2 
we have: 


ImAoj = \Aojf + \MM (24) 

Im Mj = AojM*j + MjT} 

ImTj = \Mj\M\TM ■ 

These relations are not exactly respected by perturbation theory, but are instead satisfied only to one less order in 


the expansion than kept in constructing 

the amplitude. At the 

one-loop level one has 



Im A(^J = 





ImM^^ = 


-f 



ImT^^ = 

|Mf |2 + 



For the remaining channels with I = 

J = 1 and I 

= 2, J = 

0 the UJUJ UJUJ reaction 

is elastic and the unitarity 

condition is just 






ImAij 

= \Aij? 

1^0 

(25) 

and at the NLO perturbative level. 






ImA)) 

= \Af)? 

1^0. 

(26) 


There are in all nine independent one-loop perturbative relations, that can also be obtained by applying the Landau- 
Cutkosky cutting rules and directly checked in each of the partial waves for the three reactions, providing a very good, 
non-trivial check of our amplitudes. 

Therefore the perturbative reaction matrix 

Fu = Fj^J + Fj) + ... (27) 

fulfils 

ImFjy=Fj°jFj°J (28) 

since the Ejj^ elements are real. 
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FIG. 1: Contour to apply Cauchy’s theorem in the presence of a right-hand cut (RC) due to elastic intermediate states in the 
s-channel, and a left-hand cut (LC) due to angular integration over t, ^/-channel exchanges. In the massless limit M ^ 0, the 
contour encloses only the upper half plane. The radius of the large circle is A^. 


III. THE INVERSE AMPLITUDE METHOD FOR MASSLESS PARTICLES 


A. Derivation for one channel 


The Inverse Amplitude Method (lAM) [2T] was developed for ordinary ChPT for mesons [22l [34] and it was also 
applied to the unitarization of the one-loop WBGB scattering amplitudes, at the time without a light Higgs resonance 
(see m and third reference in El). Its standard derivation is valid for one or several channels of particle pairs all 
of which have equal mass. For different masses there are technical complications (such as overlapping left and right 
cuts) that have hindered a rigorous derivation. 

In the context where we wish to apply it, for energies E ^ both masses can be taken as equal and 

negligible. Yet for massless particles, the standard derivation is also problematic, since the dispersion relation is 
thrice subtracted and the factors 1/s^ cause infrared divergences. 

Since it would be nice to have a derivation valid for massless particles, we now address a twice-subtracted dispersion 
relation that avoids infrared problems. The price to pay is that, with chiral amplitudes, the large circle at infinity 
to close the contour in the complex plane will give a contribution that needs to be calculated. As we will see in this 
section, this is feasible for elastic scattering of massless particles. 

We start by writing a twice-subtracted dispersion relation (DR) for a generic elastic partial wave amplitude A{s) 
(we suppress the I and J indices), that has both left- and right-hand cuts as shown in figure 


ds' lmA{s') s‘^ ds' lmA{s') 

TT Jo s'‘^{s' — s — ie) TT J-oo <s'^(s' — s — ie) 


(29) 


An introduction to dispersion relations can be found on refs. [ 2 S 1 E 2 . To sum up, the derivation of Eq. p^is based 


on the Cauchy theorem and on the analyticity of A{s) for Ims > 0 (first Riemann sheet), as well as on the analytic 
properties of A{s). Note that, according to Eq.[^ our computations have a left-cut (i.e., they are not analytic on the 
real axis for 5 < 0 ). So, forward dispersion relations, like those commonly introduced on text books, cannot be used. 

Because A{s) describes the scattering of Goldstone bosons, there are two simplifying properties. The first is that 
there is an Adler zero. In the massless limit this is located at 5 = 0 and guarantees A(0) = 0. Accordingly, we set 
the first subtraction constant to zero and the first term is linear in s. The second is that there are no (subthreshold, 
bound-state) poles of A{s) in the first or physical Riemann sheet (which does not make sense for Goldstone bosons 
that interact with weak strength at low energies). So only the two cuts contribute as written since A{s) is analytic in 
the rest of the upper half plane. 

We will obtain a second dispersion relation for the partial-wave amplitude expanded to NLO in the EChL, that is, 
truncated up to order A^^^( 5 ) = A^^^( 5 ) + A^^^( 5 ), that has the generic form: 


A^^\s) = Ks 

A^^^(s) = ^H(/i) + Dlog ^ + Elog . 


( 30 ) 









To derive the dispersion relation, we will first introduce the auxiliary function 


/(s) = 




(31) 


= + D log-^ + E log 

Therefore f{s) is analytic in the whole complex plane except for the left (LC) and right cuts (RC) along the negative 
and positive real axis respectively. Cauchy’s theorem provides an unsubtracted dispersion relation for f{s)\ 


1 

/(«) = - / 
^ Jo 


ds' Im/(s') 


s' — s — le 




ds' Im/(s') 
TT J_a2 s' — s — ie 


1 f ds'f{s') 
27ri Jca 


(32) 


where Ca is a circumference of radius h? oriented anticlockwise and A is an UV regulator which will be sent to infinity 
at the end (see figur e . 

Returning to Eq. ([31]), we see that this dispersion relation can easily be turned into one for 


iNLO 


5^ 

{s) = Ks + - / 
Jo 


ds'lmA^^\s') s‘^ 


s'^{s' 


/ 


ds'Im 5^ f ds'A^^\s') 


ie) TT y_A 2 s'‘^{s' — s — ie) 27 ri Jcj^ s'‘^{s' — s) 


L 


(33) 


Comparing this dispersion relation for the NLO amplitude with that for the exact amplitude A{s) in Eq. (29), we 
notice that the difference is the contribution of the circle at infinity, a term due to the divergent UV behavior of 
(X Taking now ^ s beyond the region where the amplitude is considered, the three integrals may 
easily be computed. 


s" 

TT JO 


I 

f 

J-A 


^ ds'Im 

s' 2 (s'_s_ie) = ^ a 2 

(is'Im74^^)(s') 9 „. s 


TT y_A 2 S' 2 (s' — S 

-2 r ds'A^^\s') 

Ic^ s'2(s'-s) 


i.) = 


s 

2'Ki 




r.. A 

-D log ^ 

IJ? 


r.. A 
■ E log ^ 
iB 


so that the dispersion relation for A'^^°(s) in Eq. (33) reproduces Eq. (30) 


ANLO(s) =Ks+[ B{fi) + Dlog^+Elog^ 


—s 






(34) 


(35) 


This is a consistency check of the dispersion relation and also shows its nice interplay with renormalized chiral 
couplings; the integral over the large circle trades the UV-cutoff scale A for the arbitrary renormalization scale jj. 


So far we have an elastic, exact, but not too useful, dispersion relation for A{s) in Eq. (29) and another in Eq. (33) 
for A^^^(s) which is known anyway from chiral perturbation theory. The practical use of the technique comes from 
its application to the following auxiliary function. 


w{s) 


Ms) 


(36) 


This construction has the same analytic structure than A{s) up to possible poles coming from zeroes of A{s)^ excluding 
the Adler zero (cancelled by the numerator). In addition re(0) = 0, w{s) = Ks 0{s‘^) and on the RC one has 
lmw{s) = — [A^^^(s)]^. Therefore, neglecting the possible pole contribution^, the twice-subtracted dispersion relation 
for this function reads 


s‘^ ds' lmw{s') s‘^ ds'lmw{s') s‘^ f ds'w{s') 

^ ^ ^ ^ TT Jq s'‘^ {s' — s — ie) ^ TT y_A2 s'‘^{s' — s — ie) ^ 27 ri s'‘^{s' — s) 


(37) 


^ A more careful treatment in the massive case that includes subthreshold poles found that their effect is very small, at the permille level 
or less in the physical zone m- 
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The careful choice of definition for w{s) in Eq. (36) makes possible to compute the elastic-RC integral exactly since 
Im w{s) = —K‘^s^ = Etis^ there. This is important because it is the nearest complex-plane feature to the physical 
zone (which is the upper lip of this cut, in the first Riemann sheet). 

The LC integral cannot be obtained exactly, so we choose to compute it in perturbation theory: its contribution on 
the physical zone is down by |s' — sj with s' ^ u respect to the RC, so it is small when perturbation theory deteriorates 
at 14 ^ 0. Then, it is a fair approximation to take 


lmw{s) ^ — Im 


Then one finds: 


( 5 ) c:±Ks- Ds^ log ^ - Es^ log ^ 


ds'w{s') 

' 27ri — s) 


-J. 


(38) 


(39) 


It is easy to check that this approximate integral equation is solved by w{s) = This is quite 

remarkable since w{s) in Eq. (36) is defined from the exact amplitude. Again, the only used approximations are the 
absence of poles in the inverse amplitude and the perturbative treatment of the LC integral. It stands out that, from 
the very definition of w{s), we can write down the lAM amplitude as: 




[^W(s)]^ 


yl(0)(s) - A(i)(s) 


(40) 


This I AM amplitude obtained from the ChPT expansion has many interesting properties. Eirst it has the proper 
analytic structure which, in particular, makes poles on the second Riemann sheet possible (that can be understood 
as dynamically generated resonances). Second, it is /i-invariant, depending only on the renormalized chiral constants 
a 4 , as, e, d or g that encode higher energy dynamics. It satisfies elastic unitarity, so that on the RC, 


Im A 


lAM 


_ ^lAM 


(41) 


Einally, if expanded at low energy, it coincides with the NLO-ChPT amplitude, 

^^^■^(s) = A'^l°(s) + 0(s3) (42) 

It is important to stress once more that the lAM amplitude has been obtained here by using a twice-subtracted 
dispersion relation, whereas previous derivations used a thrice-subtracted DR. Therefore we needed to carefully take 
into account the contribution of the circumference at infinity Coo, which is not present with three subtractions. 

This was necessary to avoid the infrared problems that would otherwise appear in the derivation of the lAM 
amplitude for massless particles, having all the LC and RC thresholds located at s = 0. We restate that the only 
approximations used was taking lmw{s) — ImA^^^(s) on the LC integral and assuming that w{s) has no poles, 
whereas the numerically more important RC integral is computed exactly. A posteriori these assumptions have been 
validated in low energy meson-meson scattering where the I AM method has proven to be extremely successful, as 
with a very small set of parameters it describes many different channels including their first resonances. 


B. Coupled-channel Inverse Amplitude Method 


The lAM method can also be extended to the coupled-channel case provided the masses of the particle appearing 
in the different channels are all the same, to avoid overlapping left and right cuts. This is the case here since we are 
considering the WBGB and the h particle massless. The demonstration is an immediate extension of the single-channel 
case and we relegate it to appendix [B| 

The multichannel matrix with adequate properties can be constructed from the perturbative expansion 

Fij = + ... (43) 

F(°)(s) = Ks 

F^^\s) = (B{^x) + D\og^+E\ogE\s‘^ (44) 


where now A, 5(/i), D and E have to be considered as (two by two) matrices. Eor example Kn = A, K 12 = K 21 = K' 
and K 22 = K" (Notice that K refers in different formulae to Kn or to the matrix K). Einally, Ejj is found to be 


rIAM 

^IJ 


T^i^) f 1 77 ' 


(1)^ — 1 Z7'(0) 


IJ 


(45) 
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that satisfies exact elastic unitarity on the RC 


(46) 


The various amplitudes (matrix elements of Fj^^) enjoy all the already mentioned desirable properties of the elastic 
lAM method. The coupled-channel lAM method is particular useful in the isoscalar channels (/ = 0 and J = 0 , 2 ,...) 
where the ujuj and hh channels can be strongly coupled. We dedicate section [VlT| to a detailed numerical analysis of 
the method based on Eq. (40) and (45). 


IV. THE N/D METHOD 


A. Elastic u?Cc? scattering 


The lAM is a reliable unitarization method, but to assess the systematic error introduced by approximating the 
left cut in perturbation theory, it is recommendable to compare with a different unitarization method applicable to 
the one-loop results for the ujuj and hh scattering amplitudes. A well-known alternative that we consider here is the 
N/D method. This can be applied in many different ways depending on the problem at hand. When the ujuj is purely 
elastic (J 7 ^ 0 ) the starting point is an ansatz for the scattering partial waves, from which the method is named. 


^(s) 


N{s) 

D{s) 


(47) 


where the numerator function N{s) has only a LC and the denominator function D{s) only a RC, so that A{s) has 
the expected analytical structure. Therefore Im N{s) = 0 on the RC and Im D{s) = 0 on the LC. In addition, elastic 
unitarity, ImA(s) = \A{s)\‘^ requires lmD{s) = —N(s) on the RC and we also have lmN{s) = D{s)lmA{s) on the 
LC. It is then possible in principle to write two coupled dispersion relations for N{s) and D{s). The normalization 
D{0) = 1 may be chosen by making ^"(0) = A(0), so 


D{s) 

N{s) 


1 - - 
TT 


ds'Njs') 
Jo s'{s'-s-i 
s 7° ds'D{s')lmA{s 
TT s'(s' - s-ie) 


(48) 

(49) 


More generally, one needs an n-times subtracted DR, which is useful to input the particular low-energy dynamics to 
be described: 


D{s) = 1 + hiS + h 2 s'^ + ... + hn-is^ [ ,A, ^ ^ (50) 

TT Jq s'" (s' — s — ze) 

The coupled equations for N{s) and D{s) can be solved in principle by using some recursive method. For example, 
starting from some approximate Nq{s) function featuring a LC (typically a tree level result) we can obtain Dq{s) 
by integration on the RC. Then a first approximation for the partial wave would be ^ 0 ( 5 ) = Nq{s)/Do{s). To 
continue the procedure one can now insert Do{s) in the second coupled equation to get the new a Ni{s) yielding 
^ 1 ( 5 ) = Ni{s)/Di{s) and so on. Presumably in this way it should be possible to approach as much as needed the real 
solution for some given subtraction constants, provided the original guess for No{s) is appropriate enough. Even more, 
in many cases the simplest and crude approximation A{s) No{s)/D q{s) could be considered a sensible estimate of 
the exact solution. For example, taken No{s) = A^^\s) and regularizing the integrals with IR and UV cutoffs and 
A^, one gets 


"»(•) = l + 


so that 


A(s):.. 


A(°\s) 


1 + 

We do not find this approximation satisfactory though, at least when compared with the Inverse Amplitude Method 
in section III In particular, because the equation for N has not been iterated yet, the amplitude only has a RC but 


(51) 


(52) 
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not a LC. It is unitary and depending on the UV scale A and also not compatible with the NLO result to order 

. The reason for this is that we are not yet taking into account the information coming from the NLO term 
containing the one-loop effects and the chiral couplings. However, introducing these NLO effects in the N/D method 
is far from trivial for various reasons. 

For one, it is not obvious how to choose the starting function Nq{s): remember that the NLO partial-waves have the 
general form A{s) = ( 5 ) + (s) +..., with the general form given in Eq. (15). Thus A^^^ {s) contains a logarithm 

with a LC and another one with a RC that, taken independently, are scale-dependent: the scale-independence of A{s) 
is achieved with the compensating dependence of H(/i) in Eq. (16). Thus a naive choice for No{s) featuring a LC will 
not be in general /i-invariant and that makes the N/D method less attractive. 

To solve this problem we split A^^\s) in two pieces, one having only a RC and the other only a LC and both 
/i-independent, by adequately splitting the function 


^l(s) = f log A j D52 


D + E 




Ar{s) = (AA Es'^. 


D + E 




(53) 


The cut structure is obviously as advertised, = Al{s) + Aji{s) is also trivially verified, and the scale- 

independence follows from Eq. (16). In addition, on the RC (the physical region), perturbative unitarity reads 
ImA^^^ = IirAr = The split in Eq. (|53|) is not usable in the JJ = 11 channel in the particular parameter 

case o? = b because of a coincidence^ in Eq. ( |AT8[ ), that yields E = —D. In all other circumstances the denominator 
is finite and does not give any problem. 

It is illustrative to express Al and Aji in terms of an auxiliary “loop” function 


If -s 


(54) 


This function, as the notation suggests, is /i-independent (as is easily checked). Eurthermore, it is analytical on the 
whole complex plane but for a RC. On this RC (i.e. for s = ie) we have Im^(s) = —1. Then, 


Al{s) = 7Tg{-s)Ds‘^ 
^i?(s) = 7Tg{s)Es‘^ 


(55) 


SO that perturbatively. 


A(s) = ^(0)(s) + AUs) - [^( 0 )(s)] 25 (s) + 0{s^) . 


(56) 


We have now the ingredients to apply the N/D method to the NLO computation: the useful starting point is the 
function 


No{s)=A^°'>{s) + Al{s) . 


(57) 


Notice that this function contains the LC, information about the chiral parameters and additionally it is fi independent. 

The inconvenience now is that the UV behavior of the integral for Do{s) in Eq. (48) is even worse than with the 
tree-level ansatz, since a term is included in Nq. To obtain a UV-finite integral three subtractions are required, at 
the prize of a chiral coupling of order (see appendix]^. Thus we can write: 


Dq[s) = 1 + his + h2S^ -/ 

^ Jo 


ds'[A^^\s)^AL{s)] 


/o s'^(s' — s — ie) 

As further shown in appendix 0 the N/D partial wave in this approximation can be written as 

No{s) 


A(s) - ^ 


Do{s)' 


(58) 


(59) 


It is known that in this elastic vector-isovector amplitude the NLO amplitude on the physical cut is a polynomial due to canceling 
logarithms, so the combination of chiral constants ( 0,4 — 20 , 5 ) in Eq. (A18) is /r-invariant by itself. 
















By using the Al{s) and Aji{s) definitions in Eq. (53) this denominator can also be written as: 

Al{—s)Ar{s) 


n / \ 1 ^r{^) I 1 \ /I / \ 1 ^r{s) 


2{A(0)y 
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( 60 ) 


(61) 


This amplitude in Eqs. (57), ( [5^ and ( |M| ) has many interesting properties. First it is UV finite, the IR divergences 
have been removed and it is /x independent. Second, it has the right analytical structure and it satisfies elastic unitarity 
exactly; 

Im7lN/D(s) = |^N/D(g)|2 (g2) 

on the RC. Finally it is compatible with the NLO computation up to order since: 

AN/D(s) = ^(0)(s) + ^(i)(s) + 0(3^). (63) 


All these properties are shared with the Inverse Amplitude Method. In Eq. (94) below we show that this amplitude 
converges to the lAM amplitude whenever Ar 


B. Coupled — hh channels 

Just as for the lAM, it is possible to generalize the N/D method to the multichannel case needed for the / = 0 
(J = 0, 2) cases where the ujuj state couples to the hh channel. Eollowing [39] we introduce two matrices, a numerator 
one N and a denominator D so that 


F{s) = [D{s)]-^N{s) . 


(64) 


To generalize our previous result for the single channel case, we start again from the perturbative expansion at NLO, 
Eq. (44), Again the /i evolution of is given by Eq. (16), now a matrix equation . 

Thus we can introduce the /i-independent matrix 


G{s) = -(B{fi){D + E)-^+log ^ 

TT V 1 ^ 

and the (also /^-invariant) left and right matrices 

FUs) = (^B{,,)iD + E)-^+log^'^ 
= 'kG{—s)Ds^ 

Fn{s) = {b{^){D + E)-^+\og^'^ 


Es"^ 


= nG{s)Es‘^ . 

On the RC cut perturbative unitarity reads: 


lmF^^\s) = ImFR(s) = F^°\sf = 


which implies: 



and therefore 


Fnis) = -Gis)[F(<^\s)f. 


(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 


Now we can follow essentially the same steps that we took in the single-channel case in subsection |IV A[ taking into 
account the matrix character of the different amplitudes and of K, D and E. Like in the case of the lAM, this 













13 


produces a sensible result because all particles involved, the WBGB and the Higgs-like particle, are massless and 
therefore we are not overlapping the LC and the RC in any unitarized partial wave. Then we get 

FN/D(s) = [Do{s)]-^No{s) (71) 

where 

No{s) = F^°\s) + Fl{s) (72) 

and 

Do{s) = 1 - Fr{s)[F(°\s)]-^ + |[G(s)] 277 s 2 ( 73 ) 


that can also be written as: 

Do{s) = 1 - F«(s)[F(o)(s)]-i + ^G{s)FU-s) = 1 - Fn{s)[F^°Hs)]-^ - ^Fnis)[F^°Hs)]-^FUs) (74) 
It is not difficult to check that these partial waves in Eq. 0 fulfill exact elastic unitarity on the RC, 

(75) 

and also reproduce the low-energy expansion to NLO: 

FN/°(s) =F(°)(s) + F(i)(s) + ... (76) 

Thus the F'^/’^(s) partial-wave amplitudes have all the required properties including unitarity and analyticity. They 
have a LC and RC and they can be extended to the second Riemann sheet, and in some cases have poles there that 
could be understood as resonances. 

Interesting cases where the N/D method has the advantage are those in which K = E = 0 such as the IJ = 02, 22 
waves. The vanishing of the leading-order term proportional to K makes the lAM yield zero at this order, and one 
needs the NNLO lAM or an approximation thereof, which we have not developed here but see m- However the N/D 
method can be safely applied to these situations too, as it is very easy to check since g{s) or G{s) are well defined 
even for K = E = ^. 


V. OTHER UNITARIZATION METHODS, A COMPARISON AMONG THEM, AND THEIR 

RESONANCES 

A. The K-matrix and the Improved K-matrix 

Finally we will briefly comment on some other unitarization methods which have also been considered for the 
scattering of the would-be CB in the context of the SBS of the SM. One of the most popular unitarization procedures 
is the so called K-matrix method [53] (see also [5S] for a recent review in the context of this work). The K-matrix is 
defined in terms of the S matrix as: 


I - iKl2 
I + i A/2' 


(77) 


With this parametrization S is unitary if and only if K is Hermitian. Eq. [^can be inverted to give K in terms of S': 

i(S-l) 


^ l + (S-l)/ 2 ' 
In practice the S matrix is obtained in the form of some expansion: 

s = l + S(1) s( 2 ) -h , 


(78) 


(79) 


However the truncation of this series usually produces an approximate S matrix which is not unitary. However, if we 
truncate instead an expansion of A, 




( 80 ) 




and introduce this (truncated) series into Eq. to find a new series for S, 

5 = 1 + 5 ( 1 ) + 5 ( 2 ) + ... , 
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(81) 


this is exactly unitary at any order. 

In terms of a partial-wave amplitude for some unspecified elastic process A{s)^ this amounts to the following. One 
starts from some approximate estimation Ao(s) real in the physical region and therefore not unitary. Then one defines 
the K-matrix unitarized partial wave: 


A^{s) 


^ 0 ( 5 ) 

1 - iAo(s) 


Clearly, unitarity is satisfied again in the physical region, 




1 + Ay 


(82) 


(83) 


However it is very important to stress that this K-matrix partial wave is not analytical (in the first Riemann sheet) 
and consequently it is not a proper partial wave A{s) compatible with microcausality. For example, even if Ao{s) has 
a LC, the corresponding A^{s) does not show any RC and then it cannot define a second Riemann sheet. Most of 
the criticisms on the unitarization methods appearing in the literature are based in the fact that some of them (for 
example the lAM or the N/D) can produce resonances (poles in the second Riemann sheet) but others (typically the 
K-matrix) cannot. However it is obvious, from the discussion above, that discrepancy is coming from the limitations 
(lack of proper analytic structure) of the K-matrix. The Aq (s) partial-wave is defined only in the physical region and 
cannot be extended to the whole complex plane. 

So we insist here that this naive K-matrix has no RC, cannot be extended to the second Riemann sheet, and therefore 
it cannot produce poles that could be understood as resonances. However, from our experience in unitarization 
methods in hadron physics, we know that such poles frequently appear and describe well known hadronic resonances 
[22l[34]. The original K-matrix method cannot reproduce these hadronic resonances and should be considered as less 
appropriate than other methods that are, not only unitary, but also analytical, as it is the case of the lAM or N/D 
methods. 

Nevertheless the K-matrix method can be improved as follows: we can introduce the analytical function 


g{s) = - (c + log-^ 
TT \ 


(84) 


where C is an arbitrary constant and /i is also an arbitrary scale. One interesting possibility is to define C as in 
Eq. (54) so that g{s) becomes /i independent (which is the one we will be using in the rest of the paper). In any case 
this function is analytical in the whole complex plane but for a RC. In the physical region on this RC we have: 


gis) 


- C + log^ 


■ ITT 


(85) 


and thus its imaginary part is simply 


Ims((s) = -1. (86) 

Therefore it is tempting to perform the formal substitution: —i g{s) in the K-matrix method to get what we will 
call “improved K-matrix” (IK) amplitude: 




Ao{s) 

1 + g{s)Ao{s) 


(87) 


This new amplitude is, not only unitary, but also analytical on the whole complex plane but for a RC that allows for 
analytical continuation to the second Riemann sheet, making possible the existence of poles as in the I AM or N/D 
methods. To apply this improved K-matrix method to our ujuj amplitudes, we can start by taking ^ 0 ( 5 ) = A^^\s) to 
get 


A(Q)(8) 

1 + g(s)A(^^ (s) 


A^^(s) = 


( 88 ) 
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Interestingly enough this amplitude may also be obtained from the twice-subtracted N/D method by setting in Eq. (50) 

B(iAK 


A more accurate result can be obtained by defining Ao(s) = + Ai(s) which leads to 

+ Al{s) 


or: 


A“^(s) = 




1 + 5 r(s)[A( 0 )(s) + Al(s)] 
A(0)(s) + Ai(s) 


1 - 


^h(^) _ Al{s)Ar{s) 


(89) 


(90) 


(91) 


A(0) 


(A(o')y 


This amplitude has the proper analytical behavior, is unitary and reproduces the NLO result up to order since 
Ar{s) = -5(s)(A(0))2. 

In addition, this improved K-matrix method can also be extended to the coupled-channel case, simply taking 


where again: 


F^^{s) = {l + GNor^No 

No{s)=F^°\s)+Fl{s). 


(92) 

(93) 


and G is defined in Eq. (65). 


B. The large N method 

Einally another interesting way to improve the unitarity behavior of the amplitudes is the so called large-A/" limit. 
It is based on the observation that our coset space for the electroweak SBS is SU{2)l x SU{2)r/SU{2)l+r = 
5'0(4)/5'0(3) = S^. This suggests considering a generalization to SO{N -|- l)/SO{N) = and computing the 
WBGB scattering amplitudes in the non-perturbative large- A/" limit. These amplitudes were studied in detail in na 
for the case of the minimal SM and one of their main properties is their unitarity up to NLO corrections in the 1/N 
expansion. 

However there is a limitation to the 1/N expansion as applied as an unitarization method: all channels happen to 
be I/A/'-suppressed with respect to the IJ = 00. Therefore this approximation is not appropriate to describe models 
in which other channels could be relevant, for example, those showing vector-meson dominance (such as Composite 
Higgs Bosons with low-energy W' and Z' resonances). Thus we will not consider this approach here, but we have 
shown an example of its use in |I8| . 


C. Summary of the various unitarization methods 


It has now become clear that of the several unitarization methods considered above, three stand out as acceptable, 
the version of the N/D method obtained here in sec. IV and the IK method from subsec. V A| 


the lAM in sec. Ill 


Let us gather their expressions for the elastic channels, writing them all in terms of A^i from Eq. (53) and 

g{s) from Eq. (54), for easy comparison; 


A^^'^(s) 






[A(0)(s)]2 
A(o)(s) - A(i)(s) 

A(o)(s) + A^(s) 

+9{s)Al{s) 

AW(s) + A^(s) 

^ “ AW{I) + 5fi'(s)^L(“s) 

^(°Hg) + -4L(g) 


(94) 
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All three amplitudes are IR and UV finite, /r independent, unitary, have the proper analytical structure, can be 
generalized to the coupled-channel case [see the corresponding formulae in Eqs. ( [4^ , 0 and ( [^ ] and they reproduce 
the NLO predictions of EWChPT. This attribute means that they differ from each other only at O(s^), 


= A^(s) + + O(s^) = + O(s^) = A^^(s) + O(s^) 


(95) 


Thus these three unitarization methods each provide a consistent UV completion of the low energy chiral amplitudes. 
Unfortunately, as energy grows their predictions will start differing. Then, which of them is a better description of 
reality? In principle all of them are consistent but their domain of applicability will be different. 

Eirst notice that the lAM method is the only one that does not really require the splitting of A^^^ into Al and Aji 
(or the use of the g{s) or G{s) function). This splitting is in fact in some way arbitrary, since we can always add and 
subtract a quadratic term Cs‘^ to Al and Aji respectively without changing their fundamental properties. Notice also 
that the splitting is not possible at all whenever D E = 0 (as in the / = J = 1 channel for the particular parameter 
choice a = b) and the N/D and IK methods cannot be constructed for that case. Hence, for the vector-isovector 
channel, the lAM is most appropriate. Since for D ^ E small, Al Ar^ the three methods are not expected to be 
equivalent, and we see that there are sound theoretical reasons to choose the I AM over the other two. 

Conversely the lAM method cannot be applied in the cases where K = E = 0 which happens in the J = 2 channels 
(because they start at NLO in the effective theory, so K = 0, and then perturbative unitarity forces E = K‘^ that 
also vanishes). In that case the lAM is not usable and the N/D method comes to the fore. 

In section [Vl| we will provide numerical results for the various situations to illustrate how the three unitarization 
methods work in the different channels and to try a comparison when all are applicable. Eor a brief summary, see 
table [B 


IJ 

00 

02 

II 

20 

22 

Method 


N/D 

IK 



N/D 

IK 

of 

choice 

Any 

lAM 

Any 


TABLE I: Unitarization methods usable in each IJ channel. See section IVTl 


D. Resonances 

As already mentioned, one of the more interesting properties of the lAM, N/D and IK partial waves is the possibility 
of finding poles in the second Riemann sheet under the real axis. This interest arises because these poles have the 
we can use the simple of dynamical resonances, at least when they lie close enough to the real axis in the complex s 
plane. 

Eor the amplitudes considered here the non trivial analytical behavior is coming exclusively from the logs which are 
defined in the first Riemann sheet as usual {log{z) = log(| 2 :|) -h iarg(^) with the arg( 2 ;) cut lying along the negative 
real axis). To find a pole in the second Riemann sheet, an option is to extend all the logarithms to it, through the 
simple equation 


log^^(- 2 ) = logd^l) + i[aig{z) - tt] (96) 

and then find zeroes of the denominators of the amplitudes A^^ or E^^ for coupled channels. This is the strategy that 
we followed in m- 

An alternative is to observe that given some analytical elastic amplitude A{s) = A^ {s) representing the physical 
(first) Riemann sheet, the second Riemann sheet in the quadrant under the physical region can be obtained as (see 
for example ED): 


A^\s) 


^{s) 

1 - 2iA{s) ■ 


(97) 


Therefore resonances under the real, physical s axis (the right cut) are located at points sr solving the resonance 
equation 


A{sr) + 2 ^ ® 


( 98 ) 
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so that the extension of the logarithms is unnecessary. 

The mass M and width T > 0 of the resonance can be extracted from its position, sr = — iVM. Equivalently 

we have sr = \sR\e~^^ with ^ > 0 and tan^ = q = T/M. The resonance equation (98) obviously takes a different 
form for each of the unitarization methods, which we now show in turn. For the lAM method, 


Ai^){sn) - A^^\sn) - = 0 


(99) 


whilst for the N/D method we find 

~ ^r{^r) + -^9{sr)A^^\sr)Al{—sr) — + Al{sr)] = 0 (100) 

and for the IK method, 

“ ^r{sr) + 9{sr)A^^\sr)Al{sr) — + Al{sr)] = 0. (101) 

These resonance equations are all fi independent through cancellation of their explicit and implicit (through the 
renormalized chiral parameters) dependence on /i. As expected they are different, but decreasingly so in the limit 
Al{sr) < 1, since = Ar{sr) + Al{sr). 

If we find a solution sr for some given channel IJ and some given unitarization method X = lAM, N/D, IK in 
the appropriate region M, T > 0 this solution will be a /i invariant function of the a, b and the renormalized chiral 
parameters, i.e. 


M = M(a,6,a4(/u),a4(/i)V(Ai),e(/u),5(/u);/i) (102) 


These functions trivially fulfill the observable renormalization group equations 

dM dM dM da/^ dM da^ ^ 

d/i dji da4 d/i da^ djj. 

dr dT dT da4 dT da^ 

— =-^- - H-^ + ... = 0. 

d/i djji da^ djji da^ d/i 

If we set a scale and fix the chiral couplings at that scale /io, so that 
position becomes a function of the chiral couplings evaluated at this scale only. 


(103) 


as (/io),•••, the resonance 


M = M(a, 6, a 4 , a 4 , d, e, ^) (104) 

r = T{a,b,a4^,a^,d,e,g) . 

When there is channel coupling, the amplitude matrix elements Fij (s) correspond to different reactions having the 
same quantum numbers IJ. Obviously if there is a resonance at some point sr in any of them, it should appear also 
at the same point in the rest of the matrix elements. In other words, the Fij{s) are all different as analytical functions 
but all of them have the same resonances at the same points since physically these resonances can be produced in any 
of the j ^ i reactions. 

This property is guaranteed for the three unitarization methods now at hand. This is because in all of them we 
need to invert some matrix. Thus the unitarized amplitudes Fij{s) for some given / and J contain always a common 
denominator which is a determinant depending on the unitarization method. The roots of this determinant in the 
second Riemann sheet will define the pole position for all the different processes simultaneously. 

Once we have obtained the unitarized amplitude Fij{s) by using some coupled unitarization method, and extended 
it to the corresponding second Riemann sheet (s), we can find the position of any pole (resonance) in the quadrant 
below the physical region. In the next sections we will study numerically the different channels as a function of the a 
and b parameters and the renormalized chiral couplings for the three unitarization methods considered here and we 
will compare the results obtained. 


E. Spurious resonances 

In addition to the bona-fide resonances in the second Riemann sheet, for certain sets of parameters a given unita¬ 
rization method can yield a pole in the complex s plane that lies on the first Riemann sheet. As recalled below in 
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FIG. 2: Scalar-isoscalar amplitudes (from left to right, elastic cccc, elastic hh, and cross-channel lulu hh), for a = 0.88, 6 = 3, 
and all NLO parameters set to 0 at a scale = 3 TeV. Note that, as explained on sec. |VI A| the old K-matrix method gives 
different results because its complex-s plane analytic structure is not the correct one. It will be discarded from now on. 


appendix because of Schwarz’s reflection principle, these poles always come in pairs one above and one below the 
real 5 -axis (see fig. 27 that represents the situation very graphically). 

But causality demands that the scattering amplitude be analytic in the upper half-plane, whence the pole in the 
first Riemann sheet is tachyonic. So, poles appearing in the first or physical Riemann sheet are not acceptable and 
therefore they set limits on the applicability of the method or even on the validity of the given parameter set (it might 
be that no underlying theory is compatible with the effective theory with certain parameter values). To search for 
these poles on the first sheet we just need to find zeroes of the denominators in the representations of 
and in Eq. (94). 

When we find such situation we conclude that the unitarized amplitude with the given parameters is in violation of 
causality; either there is no underlying theory that can provide such set of parameters ^ or the unitarization method 
is at its limit of validity for such set, and one can take the real part of the corresponding s where such pole appears 
as a point beyond which the theory is not applicable at all. 

Sometimes one can detect this breakup of causality in repulsive phase shifts (such as the isotensor channel) that 
vary quickly and break Wigner’s bound. 

In practice, we will consider the regions of parameter space where this phenomenon occurs as excluded. Some 
examples can be found in figures 2^ I^, ^ and ^ below. The regions where we find a pole in the first Riemann 
sheet, are automatically excluded from our parameter space. 


VI. NUMERIC COMPARISON OF THE THREE METHODS 


A. The / = J = 0 channel 


The scalar-isoscalar channel is a coupled-channel problem with the lulu and hh elastic and crossed reaction forming 
a symmetric two by two matrix. We represent the two diagonal and the off-diagonal matrix elements as functions of 
s in figure for four different methods, all of which satisfy exact unitarity. 

The three methods with the correct analyticity properties (lAM, N/D and IK) agree in predicting a scalar resonance 
that is visible in all three amplitudes between 0.8 and 0.9 TeV. The old K-matrix method gives somewhat different 
results, as known from the literature, but its complex-s plane analytic structure is not the correct one, as visible in 
Eq. ( [ 8 ^ . We therefore discard the old K-matrix method from now on. 

The other three methods are practically in perfect agreement up to the first elastic resonance and they start deviating 
quantitatively only f or hig her energies. The reason that there is good agreement between the various methods was 
discussed under Eq. (101): since we have set the NLO terms to 0, Al is small, and the three resonance equations 


become dominated by the tree-level and right-cut parts of the amplitude, which suggests similar masses for all the 
methods. 

Note that in [42] we have shown that the resonance found in figure appears even if we set a = 1 (its SM value 
with one Higgs): it is sufficient that the coupled-channel dynamics is strong through — 6 7 ^ 0 for it to appear. 


5 


Other authors speak of “negative width resonances”, presumably because of Eq. ( |l09| , but we do not favor this concept, as it seems a 
linguistic contradiction in terms. 
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FIG. 3: Vector-isovector partial wave. We have taken a = 0.88 and b = 1.5, but while for the left plot all the NLO parameters 
vanish, for the right plot we have taken a 4 = 0.003, known to yield an lAM resonance from the work of the Barcelona group [6]. 
Note that the N/D and K-improved methods are not reliable in this channel, as explained below on section VIB They are 
included to show the lack of agreement with the I AM. 


Moreover, with the values chosen to prepare the figure this — 6 i s negative, so the cross-channel amplitude Mj 
shown in the right-most plot is also negative as dictated by Eq. ( |A22| ). At last, observe that the resonance appears in 
all three elastic or inelastic amplitudes in the same position (though of course, with different shapes due to different 
backgrounds). 


B. The / = J = 1 channel 


The comparison between the three methods lAM, N/D and IK for the vector-isovector channel is shown in figure 
First we set all the NLO parameters to 0 (left plot). Clearly, there is no good agreement between the three lAM, N/D 
and IK methods. Moreover, if we introduce one NLO counterterm with an appropriate value to generate a resonance 
in the lAM, here a 4 = 0.003 as an example (right plot), the N/D and IK methods do not react in the same way as 
the former, and fail to yield a vector resonance. 

In order to understand the discrepancy found in this elastic channel we notice that the possibility of defining the 
N/D and the IK methods depends dramatically on having D E 0 since otherwise we cannot define the g{s) 
function in Eq. (54) nor the Al and Ar splitting in Eq. (53). 


But here comes the coincidence, in the / = J = 1 channel we have 


D, 


-El 1 — 


(96)^7r^u^ 


(a2-6)2 


(105) 


which vanishes for = b. This is in particular the case of the SM where a = 6 = 1, which is not very important for 
our discussion because there are no strong interactions to start with. More importantly, = 6 is also satisfied by the 
Higgsless electroweak chiral perturbation theory, characterized by a = 5 = 0. This situation is already ruled out by 
the discovery of the light Higgs-like particle, but it is still interesting because it is equivalent to two-flavor low-energy 
QCD in the chiral limit with v playing the role of and the WBGB being the pions. 

Within a = 0 = 6 , we know that a vector resonance (the p) appears in the spectrum (because we can look up the 
answer in QCD), and know what the low-energy parameters are, with good approximation. Figured shows the result 
of the calculation with the lAM (solid line). We have there taken a = 6 = 0 and a 4 = — 2 a 5 = ^^ 2 , the large-A^c 
prediction for these NLO constants (others taken to 0). The p vector-isovector resonance then comes with reasonable 
parameters (to see it, substitute v = 246 GeV by / = 92 MeV in the scale; this amounts to rrip 2.1 TeV ^ 0.79 GeV, 
just slightly above the actual 0.775 GeV in the hadron spectrum). 

The other lines in figure have been computed by increasing a towards 0.88, the value taken for figure One sees 
without doubt how the QCD-like resonance becomes narrower and lighter (this depends on the interplay of a with the 
NLO parameters a 4 , as), matching the calculation of figure]^ We find that the IK and N/D methods fail to provide 
a resonance. Therefore, the lAM is the method of choice for the vector-isovector channel, given that the other two 
fail at least over the o? = b parameter election, while the lAM yields a resonance that can be continuously matched 
to the one we know is there for that parameter set. 
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FIG. 4: We show the vector-isovector resonance with NLO a 4 , as parameters taken from large-iVc QCD, b — o? and a as shown 
in the legend. The right-most solid line is the rescaled QCD case, towards the left we approach the EWSBS with a Higgs, 
where the resonance is narrow and relatively light for these a 4 , as. 


The resonance may be exactly fit to data with an adequate choice of the a 4 and chiral parameters to adjust its 
mass and width. Beyond trial and error, an elegant method is to couple the resonance to the Chiral Lagrangian in 
a chiral invariant way and then integrate the resonance at tree level as done for example in [33] (see also the early 
treatment by [H] and the more formal one in [45], as well as that in the context of Composite Higgs Models in [29]). 
The tree-level chiral couplings obtained take the general form 




VTl^-treey 


(106) 


where i = 4 , 5 , 774 = —775 = 127 r and 7^^®® = with F^^®® being the tree-level vector-resonance 

parameters. Thus the tree-level term induced by the resonance is 

A\Tis) = (107) 

where the 794 and 795 constants are obtained from Eq. ( pT| ), = Bq p 4 a 4 {fi) and are given by 

P 4 = l/(247r7;^) and p^ = — 2 ^ 4 . Following [ 43 ] we can now obtain the contribution to the renormalized chiral 
couplings induced by the resonance by matching the O(s^) tree level amplitude with the NLO result at the point 
s = 2 i. e. 


This identification leads us to 


^tree(^tree 2^ ^ ^). 




Bo 

P4+P5 


for z = 4,5. Therefore we get 


(s) = 


37* 


27Vftree 4 


Dll log 


T^tree 2 


Ell log 


—s 


]\ptree 2 


(108) 


(109) 


( 110 ) 


Then the lAM resonance equation 
resonance limit 7 = F/M ^ 1: 


(99) leads us to the 


second Riemann-sheet resonance parameters in the narrow- 


(^IAM )2 

pIAM 


i^ll 

Bll(Mt>'ee) 

^2^^IAM 

Bll(M‘"'««) 


( 111 ) 
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FIG. 5: Scalar-isotensor amplitudes for a = 0.88, 6 = and the NLO parameters set to 0. All three unitarization methods 
agree qualitatively and with the perturbative amplitude too, as loop corrections are small. Here we plot both the imaginary 
part (top set of lines) and the real par t (bot tom set). That the real part is negative reflects the repulsive interaction in this 
channel given by —(1 — a^) < 0 in Eq. (A19). 


which implies the M^^®®-independent result = Kn or: 


pIAM 




(l-a^) 


( 112 ) 


which is recognizable as a version of the KSFR relation [461147] (slightly generalized to a 7 ^ 0). This is here a restriction 
arising from the constraint of exact unitarity, that has been discussed in m and references therein and is a non-trivial 
relation between three observable quantities. 

Also we have the equation: 


(113) 

which relates the resonance parameters with the tree level ones. This is a very consistent result showing that the I AM 
method properly predicts a vector resonance whenever , F^^®® > 0 , in which case the chiral parameters receive a 
contribution and may be dominated by a vector resonance. For example implies F^^^ = (3/2)F^^®® 

which is a quite reasonable result taken into account the tree-level nature of the vector field integration performed to 
estimate the chiral parameters. 

However the N/D and IK unitarization methods fail to predict this resonance for the appropriate values of the 
chiral parameters. First they are not even defined for a = b. For a 7 ^ 6 but still in the parameter region close to the 
SM where a ^ 6 ^ 1 we have i^n + ^ 0. In this case the methods are well defined but then Al ^ Aji which 

means that the lAM method is very different from the N/D and IK methods. Thus, as the lAM method works pretty 
well in this channel according to the previous discussion, we have to conclude that the other two methods are not 
appropriate to describe the vector channel. 


M 


lAM 


M' 


tree 


/2riAM 

(3f^ 


C. Scalar-Isotensor J = 0, / = 2 channel 

We now consider the isotensor channeHwhere a resonance, if there ever was one, would distinctly appear for example 
in equal-charge spectra). Figure^ shows the resulting amplitude for a = 0.88, b = o? and all NLO parameters 

set to 0 . 

We plot both the real and the imaginary parts of the three unitarized amplitudes and obtain a moderately weak, 
repulsive partial wave that does not bind a resonance (as seen from the negative real part). All three unitarization 
methods give a consistent picture: the unitarized interaction has a slightly larger real part and slightly smaller 
imaginary part than the (unitarity-violating) perturbative one. 

In figurein turn we plot the same isotensor amplitude for a = 1.15. Now the real part has opposite sign (attractive 
interactionjand grows more rapidly, with all the unitarization methods agreeing and once more tracking perturbation 
theory until about the end of our energy interval at 3 TeV. 














22 



FIG. 6: Scalar-isotensor amplitudes for a = 1.15, b = and the NLO parameters set to 0. All three unitarization methods 
agree qualitatively once again, even though now the amplitudes are strong. The real part (corresponding to the set of lines 
larger at low-F^, since it receives a tree-level contribution unlike the imaginary part) is now positive because of the sign reversal 
of (1 — a^) respect to figure]^ 



FIG. 7: Tensor-isoscalar amplitude for a = 0.88, b as shown, and the NLO parameters set to 0. The amplitude is real. 


D. Tensor isoscalar channel with J = 2, / = 0 

In hadron physics there is a well known /2(1270) resonance that is broad and visible in tt+tt” spectra. Its mass is 
well above the 775 MeV of the vector p, which is natural because the d-wave is smaller than the p-wave due to the 
suppression factor near threshold. 

In figure we show the tensor-isoscalar channel in perturbation theory, which is indeed small, with all the NLO 
parameters set to 0, and a = 0.88, b as shown in the figure legend. This is once equal to to show the elastic 
amplitude, and once equal to a^/2 to see the other, inelastic and hh amplitudes. All are of course real and quadratic 
in s (because Kq 2 = 0, the LO 0{s) vanishes). 

Next we show, in figure the comparison between the N/D and IK method in unitarizing the partial wave with 
/ = 0, J = 2. The I AM method vanishes and cannot be used without information from NNLO, because here the LO 
in perturbation theory is zero {Kq 2 = 0). 

In the left plot we have set a = 6 = 0 and = —2a^ = as in figured The IK method clearly shows, and 

the N/D method is suggestive of, a QCD-like /2 resonance (rescaling again v = z46GeV to = 92 MeV, the 3.5 TeV 
resonance mass becomes 1.3 GeV, in very good agreement with the experimental 1.27 GeV /2 resonance in the hadron 
spectrum; and this with no free NLO parameters, since they are taken from large-Vc). 

In the right plot we have now increased a = 0.88, with b = o? still fixed to avoid the coupled-channel situation. The 
resonance is seen to become lighter and narrower, and both unitarization methods qualitatively agree in predicting 
the resonance though the mass is slightly different. 

If we now lift the b = o? requirement, because this is an isoscalar channel the hh system becomes coupled to ujuj. 
Then the resonance should be visible in both particle spectra, and also in the channel-coupling amplitude; all three 
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FIG. 8: Comparison of the two available methods of unitarization for the isoscalar-tensor channel / = 0, J = 2 with 6 = 
(only one channel). The as constants have been fixed to their values in large-A^c gauge theory, so the left plot with a=0 
reproduces the QCD situation with a broad, heavy / 2 -like resonance. The right plot for a = 0.88 shows how this becomes 
narrow. Both methods agree well. 
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FIG. 9: Isoscalar-tensor amplitudes (imaginary parts) for a — 0.88, b = a^/2, and the NLO parameters set to 0. From left to 
right: elastic lulu, elastic hh and cross-channel amplitudes. 


are shown in figure where the now inelastic resonance is clearly visible. 

Its mass is very similar to the purely elastic case, and both unitarization methods continue being in qualitative 
agreement. 

We use the opportunity to show the appearance of this resonance also as a consequence of the channel coupling 
induced by the parameter e of the effective Lagrangian. The lAM below does not capture the tensor channel, and the 
scalar one that the lAM does capture is only sensitive to the combination d + e/3 which does not allow to disent angle 
the two parameters. To see the separate effect of e we need to examine the tensor channel^ as seen in Eq. (A23), and 
this can be carried out with the N/D or IK methods. We show the result of the analysis in figure To prepare the 
figure we have taken a = 0.95 and h = a?{2. If all the NLO parameters vanish, there is no low-energy resonance in 
this tensor-isoscalar channel. Adding e at the level of3 — 4x 10“^ or more causes a resonance to enter the low-energy 
region. 


E. Tensor-isotensor channel with J = 2, / = 2 


The last partial wave that does not vanish at one-loop order in perturbation theory, and that to our knowledge has 
not been considered in the literature, is the tensor-isotensor channel. Here again K 22 = E 22 = 0 so that the amp litude 
in perturbation theory is real for physical energy. The non-vanishing constants, B 22 and D 22 are given in Eq. (A21) 
below and the amplitude is drawn in figure in perturbation theory. 


6 


This arises naturally because the d^hd^h contraction that multiplies d in Eq. (|^ is a scalar, while the d^hd^h one that accompanies e 
has both scalar and tensor components. 
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FIG. 10: The tensor-isoscalar J = 2, / = 0 coupled channels analyzed with both IK (dashed lines) and N/D (solid line) methods 
can show a resonance induced by the parameter e. 
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FIG. 11: The real tensor, isotensor I — J — 2 amplitude in NLO perturbation theory for a — 0.88 and two values of h. 

Moreover, figure [l^ shows this computation in perturbation theory for the case h = o? together with the isotensor- 
scalar one and also the two isoscalar amplitudes. Comparing those of equal I we see that larger J is suppressed below 
Attv ^ 3TeV (more so for the scalar channel, since the scalar-isoscalar amplitude is strongly interacting). Curiously, 
for J = 2 the isotensor wave is stronger than the isoscalar one. 

The unitarization of the J = I = 2 channel is not possible in the lAM method because K 22 = 0, but both IK and 
N/D methods concur in the presence of a resonance, as seen in figurewhen the NLO parameter is large enough. 
It is worth remarking that, for a given a 4 , mu < m 22 so that having this resonance in the 2-3 TeV region entails the 
presence of the vector-isovector (p-like one) in the 1-2 TeV energy interval. 

As we have established that the convergence of the partial wave expansion is very good by comparing the J = 2 and 
J = 0 amplitudes, and that the order of the spectrum of resonances is the natural one, with those of lower angular 
momentum appearing at lower energy, we concentrate in the following on the three cases that are accessible to the 
NLO-IAM, the 00, 11 and 20 channels; only the first one requires the coupled-channel treatment. 


VII. SYSTEMATIC NUMERICAL STUDY OF THE lAM 


In this section we undertake the systematic study of the lAM with the help of a computer. The calculations are very 
straightforward and involve simple algebraic formula (no integrations, as the dispersion relation has been analytically 
solved) and the inversion, at most, of dimension-two matrices. The I AM cannot handle, without NNLO information, 
the higher partial waves with J = 2 or beyond, but we have seen in figure that, under natural conditions, these 
are quite smaller in the low-energy region. For the three dominant low-energy amplitudes, the lAM based on NLO 
perturbation theory is reliab le and powerful, so we proceed with it alone. 

First, in subsection VIIA we address the one-channel I AM in Eq. (40) for the WlWl elastic scattering, with the 
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FIG. 12: Moduli of the isoscalar and isotensor NLO perturbative amplitudes theory for a = 0.88 and h = o?^ showing good 
convergence of the partial wave expansion in the low energy region (the J = 2 waves are much smaller than the two J = 0 
waves). 
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FIG. 13: Tensor-isotensor resonance as function of the NLO parameter for the IK-matrix (left plot) and the N/D method 
(right plot). 


help of the equivalence theorem, of course. This involves setting 6 = 0 and studying the behavior of the amplitudes 
upon varying each of the acti ve par ameters a, a 4 and a^. These results are just reassuring as they are known to a 
large extent. Then subsection VIIB addresses the coupled channels, by means of Eq. (45) and it is here that we make 
a totally new contribution. 

One of our findings is a coupled-channel resonance akin to the low-energy a meson but that can be generated by 
purely ww — hh interactions independently of the elastic potential strength between two ws or two hs. We have chosen 
to highlight this curious object in a companion letter [12] so we do not focuse on it so much here. 


A. Purely elastic scattering with b = o? 

The current 2cr bounds on the a parameter are, from CMS, a G (0.88 — 1.15), and a G (0.96 — 1.34) from the 
ATLAS collaboration [52]. We will take as reference a fixed value of a = 0.95 with NLO parameters set to 0, and later 
exemplify the sensitivity to each parameter (a is better chosen different from 1 because of the factor (1 — o?) that 
enters the leading order amplitudes). In any case, the sensitivity to a is displayed in figure pM} Generally speaking, 
for a < 1 (left plot) there is a broad scalar resonance akin to the cr in hadron physics, and the other channels are 
nonresonant. For a > 1 we can see a different situation in which the scalar strength significantly diminishes, but 
instead the isotensor wave becomes strong and possibly resonant (because the factor 1 — changes sign, so its 
normally repulsive amplitude becomes attractive). 

We now take the middle plot in figureand add an NLO term proportional to either a 4 or as, with the outcome 
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1=0, J=0 - 1=1, J=1 . I=2, J=0 - l=0, J=0 - l=1,J=1 . 1=2, J=0 — 1=0, J=0 - 1=1, J=1 . 1=2, J=0 





FIG. 14: Moduli of the lowest elastic ojoj ojoj partial waves in the I AM for h — (no coupled channels) as function of a. We 
will take the middle plot as reference for the parameter exploration in the next graphs. From left to right, a = 0.75, 0.95, 1.25. 


1=0, J=0 - 1=1, J=1 = l=2, J=0 - l=0, J=0 - l=1,J=1 . 1=2, J=0 




FIG. 15: Moduli of the lowest elastic uu ujuj partial waves in the I AM for h — c? (no coupled channels) showing the effect of 
a 4 (left) and as (right) both positive and alternatively equal to 0.002. Here a — 0.95. We see a light scalar-isoscalar resonance, 
a vector-isovector resonance around a TeV in the left plot (that moves to higher masses for smaller values of the positive a 4 
that induces it), and an inconspicuous isotensor amplitude. 


plotted in figure [T^ The effect of a 4 of order 10“^ (left plot) is to produce a very narrow vector-isovector resonance, 
and narrowing plus making lighter the scalar-isoscalar one. The effect of (right plot) at this same level of intensity is 
only dramatic in the scalar-isoscalar channel, while the vector one remains of moderate intensity and hardly resonant 
at all. This is in agreement with the observation in [6]. The vector resonance induced by positive a 4 can also be seen 
in the scattering phase shift in figure pT) The left plot shows the phase motion in the three lowest-F^ channels with 
all NLO parameters set to 0. No resonance is seen, in agreement with the middle plot of figure The right plot 
shows clear resonant phase motion corresponding to the resonances in figure where we study the effect of both 
a 4 and a^. The good agreement with [6] is remarkable, both works agreeing on the appearance of a pole on the first 


- 1=0, J=0 - 1=1, J=1 . 1=2, J=0 



1=0, J=0 - 1=1, J=1 


-- 1=2, J=0 


FIG. 16: Scattering phase shift of the lowest elastic ujuj ujuj partial waves in the lAM for h — — 0.95^ (no coupled 

channels), a 4 = 0 (left plot) and a 4 = 0.002 (right plot). We can see how indeed the addition of an a 4 at the level of 10“^ 
generates phase motion crossing 7r/2 in the right plot corresponding to a resonance in both the scalar and vector channels. 
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FIG. 17: From left to right, isoscalar {IJ — 00), isovector (/J = 11) and isotensor {IJ — 20) channels in elastic uu ojoj 
scattering. For a = 0.90 (different from our base scenario so we may compare with other groups), 6 = a^, we show the a 4 -a 5 
parameter map, setting the other NLO parameters to zero. Note the appearance of a pole on the first Riemann sheet for 
/ J = 20 and negative enough values of both a 4 and as. The comparison with ref. [6] is very satisfactory. 



FIG. 18: Example pole of the isoscalar elastic amplitude with b — o? (only the ujuo uju channel is active), a = 0.95, a 4 = 10“^, 
and all other NLO parameters set to 0. Pole in the second Riemann sheet (below the physical, real-s axis highlighted in bright 
yellow online). The lower (salmon online) and upper (blue online) surfaces are, respectively, the first and second Riemann 
sheets. 


Riemann sheet in the isotensor channel for negative enough values of either a 4 or a^. 


This feature is shown in figure is in full agreement with the results of [6] and, as discussed in subsection V E 


excludes this parameter space within the I AM. The computational method to produce this and the following maps in 
parameter space is described in appendix [D| 

In figure [l7| we call the experimentally disfavored regions so because poles appear with \s\ < (700GeV)^ (scalar- 
isoscalar and isotensor channels) and (1.5TeV)^ (vector-isovector channel). 

The vector-isovector channel is here exceptional in that the two variables enter with opposite signs, in the combi¬ 
nation a 4 — 2a5, see Eq. whereas in all other four NLO amplitudes they come with equal sign. Thus, the slant 

in the middle plot is opposite to the other two. 

For broad swipes of a 4 -a 5 parameter space the lAM predicts either isoscalar or isovector resonances or both. In 
figure p!8| we show an example of a pole in the second Riemann sheet of elastic ujuj scattering in / = 0, the Aqo partial 
wave for one channel only {b = o?). 

Therein the continuation to the second Riemann sheet has been obtained with Eq. (97) and the resonance appears 
as appropriate below the real, physical s-axis (bright yellow line). This pole corresponds to the scalar lAM resonance 
shown for physical s in figure (blue solid line there) though a 4 is somewhat smaller here. This serves as illustration 
of the pole structures in the complex plane (unstable particles or resonances) that accompany our resonant shapes 
for physical s. 

A lot of the a 4 -a 5 parameter space represented on [6] is experimentally disfavored because the mass-range where 
the resonances appear is being covered by LHC data [3], with none found yet, though such experimental bounds are 
not very strong because the couplings between the new resonances and the detected SM leptons are quite arbitrary 
(from the effective theory point of view), so it is difficult to interpret the bounds beyond particular models. 
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FIG. 19: From left to right, isoscalar {IJ — 00), isovector (/J = 11) and isotensor [IJ — 20) channels, o? — b vs. a 4 . Note the 
presence of poles on the first Riemann sheet for certain region of the a > 1, a 4 < 0 parameter space. 


ww->ww - ww->hh . hh->hh - ww->ww - ww->hh . hh->hh - ww->ww - ww->hh . hh->hh 





FIG. 20: Moduli of the lowest (/ = J = 0) partial waves in the I AM for h ^ = 1 (all the strong dynamics comes from the 

coupled channels). From left to right, 1, 5 = 2, 5 = 3 (the first and third are almost equal since they are symmetric 

respect to 5 = 1). A scalar resonant structure is apparent for = 1 TeV; because more extreme values of 5 lower its mass, we 
are able to give a bound on the value of 5, that must be roughly contained in (—1,3), as explained in the companion letter [42]. 
We will take the middle plot as reference for the parameter exploration in several of the following graphs. 


On fig. 19 the simultaneous effect of a (with o? = 5) and a 4 is shown, again swiping the parameter space looking 
for resonances. Note the presence of a resonance on the first Riemann sheet in the isotensor channel even for a < 1 
and sufficiently negative values of a 4 . For a > 1 (and 5 = a^), there is no resonance on the first Riemann sheet. For 
a < 1, we can find a pole in both the isoscalar and isovector channels. For a > 1, only an isotensor resonance is to be 
found. 


B. Scattering in the presence o^ h ^ o? 


Setting h a? = 1 opens the inelastic scattering ujuj hh channel in the absence of elastic strength. Figure 
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shows the dependence on 5. Almost all our computed perturbative amplitudes are symmetric around 5 = o? = 1 


(see sec. A3), with the exception of the scalar-isoscalar ujuj hh channel-mixing Mq partial wave in Eq. A22 


this asymmetry then appears in other channels due to the unitarization (a way of thinking of it is with the image 
of resumming perturbation theory), but the effect is small, so that the left and right plots are quite similar. The 
scalar-isoscalar resonance shown is very interesting and the object of focuse of the accompanying letter [42] . 

shows the lowest elastic ujuj ujuj partial waves in the presence of a 7 ^ 1 (as well as 5 7 ^ a^), so there is both 


Figs. 
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elastic and inelastic potential strength. The scalar resonance is then more similar to the standard QCD a resonance. 

A novelty is the appearance of a pole on the second Riemann sheet of the isotensor channel for a = 1.25, 5 = 1.1. 
This is very much unlike QCD, where the isotensor channel is weak and repulsive; while there is no 7 r+ 7 r+ resonance 
in the hadron spectrum, this is still allowed by current constraints on the one. 

However, as we show on figurethis case with a > 1 is quite critical, because most of the parameter space features 
an isovector pole on the first Riemann sheet, so that much of this parameter region must be ruled out or declared 
beyond our validity range. Only a small part of the a > 1 parameter space shows an isotensor pole on the second 
Riemann sheet while excluding an isovector pole on the first one, and simultaneously remains out of experimentally 
disfavored values of a. 

On the other hand, the behaviour for a < 1 is more standard, showing a resonance on the second Riemann sheet 
only in the isoscalar channel. This resonance is quite broad, and only becomes experimentally disfavored for relatively 
large values of — 5. 
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1=0, J=0 - 1=1, J=1 -------- l=2, J=0 - l=0, J=0 - l=1,J=1 . 1=2, J=0 





FIG. 21 : Moduli of the lowest elastic cjcj ^ ujuj partial waves in the I AM for 6 / 7 ^ 1 (strength from both elastic and 

coupled-channel dynamics). Left plot: a — 0.75, h — 0.9, showing much strength in the scalar channel, presumably due to a cr 
resonance. Center: a = 1.25, b — 1.1, showing a pole on the second Riemann sheet in the isotensor channel, clearly seen also 
on the right plot in the complex plane. 
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FIG. 22: From left to right, isoscalar (/J = 00, isovector (/J = 11) and isotensor {IJ — 20) channels. Note the presence 
of a pole in the first Riemann sheet of the isovector channel in quite some of the parameter space with a > 1. All the NLO 
parameters are set to zero. 


The d and e parameters are studied on figures and respectively. However, note that they appear in the 
combination d + (e/3) on the lowest partial wave {IJ = 00), so the lAM applied to any future strongly coupled 
resonance would be insufficient to separate them and one would need to resort to the J = 2, / = 0 resonance in 
figure above to obtain e independently of d. 

We concentrate now on the / = J = 0, a = l,5 = 2 case, which has an isoscalar pole on the second Riemann sheet. 
hh is shown on figures (right) and 


A peak on ujuj - 
operators ujuohh 
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This is expected, since d and e accompany four-particle 


In fig. [^we see that for positive values of d or e, the isoscalar pole weakens and then disappears. 
But for negative values, a pole on the first Riemann sheet emerges. The case of d = —0.01 shown in figure|^is curious 
because there is no pole on the first Riemann sheet below 3TeV so we should not a priori reject all that structure 
in the corresponding plots of figure including a zero of the amplitude at high energies. Of course, we should be 
cautious: perhaps, for these small negative values of d the pole simply moves to higher energies and we should not 
trust the computation (or discard negative d altogether). 

Finally, we study the dependence of all amplitudes on the g parameter (the only one that we have kept from the 
pure Higgs scattering sector, as it is needed to renormalize our amplitudes). 

It most directly produces an enhancement of hh hh scattering that starts at NLO, as can be seen in figure [26l 
since it comes from a {d^hd^h^ term in the effective Lagrangian. 

In figure we study the parameter combination a = 1 and 5 > 1.5 together with a varying g, so we see the 
interplay of the channel coupling with the Higgs-sector dynamics. 

We find a proper isoscalar pole on the second Riemann sheet for positive g. If either g or b are somewhat large, the 
isoscalar resonance enters the experimentally disfavored zone where LHC data are having an impact. 

On the contrary, negative values of g introduce a pole on the first Riemann sheet, so we must exclude those. 
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ww->ww - ww->hh . hh->hh 



ww->ww - ww->hh hh->hh 



ww->ww - ww->hh hh->hh 



ww->ww - ww->hh hh->hh 



FIG. 23: Sensitivity to d. We depict the lowest (/ = J = 0) partial wave in the I AM for 6 = 2 7 ^ = 1 . Left: moduli of the 

amplitudes with d = 0.01 (top) and d = — 0.01 (bottom). Right: real (top) and imaginary (bottom) value of that partial wave 
for d = — 0 . 01 , where we see that the channel-coupling partial wave is analytic but has a zero. 


-. ww->ww - ww->hh . hh->hh ww->ww - ww->hh hh->hh 




FIG. 24: Moduli of the lowest (/ = J = 0) partial waves in the lAM for h — 2 ^ — 1. Left plot: e = 0.01. Right plot: 

e = —0.01. The result is similar to fig. |23| because, of course, this channel depends only on the parameter combination 
d + (e/3), which serves as a check. 
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FIG. 25: Scalar-isoscalar channel {IJ — 00), with a = 1 , 5 = 2 . Left: d-e parameter map looking for poles. Right: imaginary 
part of the elastic ujuj scattering (^d — e — —0.005). The isovect or an d isotensor channels, not shown, have no poles in the 


region of interest (|s| < (3TeV)^). As discussed above in subsec. 
sheet (and a conjugate pole that is outside our circuit). 


VE 


the black region contains a pole on the first Riemann 



FIG. 26: Dependence on that we find weak for natural values thereof. Displayed are the moduli of the lowest (/ = J = 0) 
lAM partial waves for 5 = 2 7 ^ = 1. Top panel: from left to right, g — 0.002, g — 0.005, g — 0.01. Bottom panel: negative g 
values of equal magnitude. 
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ba 



I I No pole 

Exp. allowed 
Disfavored 
Pole 1st 



FIG. 27: Left: map of the h-g parameter space, seeking poles in the isocalar channel {IJ — 00), with a = 1 fixed and 
the remaining NLO parameters set to zero. (The isovector and isotensor channels have no poles in the region of interest, 
\s\ < (3TeV)^.) In the black regions there are two poles above and below the real axis on the first Riemann sheet, and we 
capture at least one with Cauchy’s theorem, excluding the corresponding parameter swath. Right: explicit plot of these two 
poles for fixed parameter values b = 2.4, g = —0.08 (plotting again the imaginary part of the elastic lulu scattering). 
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VIII. SUMMARY AND DISCUSSION 


In this article we have presented a thorough study of the unitarization of the Effective Lagrangian describing the 
Electroweak Symmetry Breaking Sector in the TeV region. 

The Effective Lagrangian in the massless limit has seven free parameters, namely a and b that respectively provide 
elastic ujuj ujuj and cross-channel ujuj hh strength at LO, and five more at NLO: the elastic and (inherited 
from the old Electroweak Chiral Lagrangian), d and e (that couple the two channels at NLO) and g (in the pure 
hh hh sector). This is the minimum number of parameters necessary to obtain a renormalized theory at NLO for 
massless uj and h bosons. The parameter set, the combinations in which they appear, and the experimental reactions 
useful to extract them are summarized for convenience in table im 

We have discussed five unitarization methods, aiming at classifying their respective strengths and weaknesses. 
We argued that three of them satisfy all desirable properties (describe several IJ channels, produce unitary and 
analytic amplitudes, are independent of the renormalization scale, and agree with perturbation theory at low energy) 
and provided explicit constructions for them based on exact (elastic) dispersion relations. These are the Inverse 
Amplitude Method, that we have studied at length, the N/D method and the Improved K-matrix method, that we 
have also assessed. All three have been compared. 

The three methods are applicable to the / = J = 0 coupled-channel partial wave, and to the exotic / = 2, J = 0 
UJUJ channel. Eor any given set of parameters in the Lagrangian, the three methods are in qualitative agreement. In 
particular, they all produce a a-like resonance when the interactions become strong, and the mass values obtained 
agree to within a few percent, which is quite remarkable and means that the model dependence is well controlled by 
imposing all the necessary theory constraints. 

We have also unveiled a coupled-channel /o-hke scalar-isoscalar resonance that appears even if a 1 as long as 
b is large enough (to provide coupled-channel strength). We have written a companion letter to this already long 
article highlighting this resonance. We only remark here that, though the LHC starts imposing relatively significant 
constraints on the a parameter, it has not made much progress of substance in constraining 6, so this coupled-channel 
resonance is one of the most interesting strongly interacting objects that can be sought for at the LHC run-H and 
beyond, because it may appear at relatively low-energies of 1 TeV or less (because of its somewhat large width). 

In the / = 1 = J channel (covering for example the W' and Z' bosons associated to Composite Higgs Models, as 
long as they be strongly coupled to ujuj) the lAM is the method of choice because the other two cannot be constructed 
in a renormalization-scale invariant way. 

Einally for the two channels with J = 2 (where in particular / 2 -like resonances might appear, as well as exotic 
ones in same-charge combinations) the lAM cannot be constructed with NLO amplitudes (because the lowest 

order is for these), but the other two methods do work and are in qualitative agreement. 

We have provided extensive numerical analysis of all these amplitudes, for physical, real values of 5 , as well as into 
the complex s plane. Therein we have searched for poles of the scattering amplitudes in their second Riemann sheet 
to be interpreted as resonances, as well as poles in the first Riemann sheet that exclude certain regions of parameter 
space. We have then drafted bidimensional maps of the parameter space showing whether the poles are likely to be 
excluded by LHC searchers, are in violation of causality, or are still viable resonances that can be searched for, and 
they agree with those in prior literature where available. 

To conclude, we believe that we have made a substantive contribution to the discussion of possible strongly- 
interacting extensions of the Standard Model, which are the currently most natural scenario with the found particles 
W, Z and h in the Electroweak Symmetry Breaking sector. Thus we have extended previous works done long before 
the discovery of the 125 GeV Higgs like boson, that did not include it, as for example those in [53]. Of course, it 
can still be that the SM exhausts TeV-scale physics, in which case the parameters of the effective Lagrangian become 
a = b = 1 (all the NLO ones vanishing). Or it can also be that some of them only slightly deviate from the SM 
values; this could be suggestive of weakly coupled resonances, as per Eq. (106) and the theory would be unitary far 
from saturation. But the strongly interacting regime remains the bulk of the parameter space to be explored by the 
LHC run-H. 
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Parameter 

Combination 

Simplest reactions 

Expt. extraction 

Resonance type 

a 

1 — a^ 

^00,^115^20 cx: (1 — a^)s 

Low-E WlWl s-wave 

cr-like 

b 

— b 

Mo oc {0? — b)s 

Low-E hh s-wave 

Coupled-channel fo 

0,4, as 

2a4 + as 

Ao2 oc [(...) + (... )(2a4 + as)]s^ 

A20 oc (...) s + [(...) + (... )(2a4 + as)]s^ 

Low-E WlWl d-wave 

Low-E WlWl d-wave 

/2-like 

Exotic VE+VE+ 

p-like 

a 4 + 2as 

A22 oc [(...) + (... )(a4 + 2as)]s^ 

7a4 + lias 

Aoo oc (...) s + [(...) + (... )(7a4 + llas)]5^ 

a4 — 2as 

All oc (...)s+ [(...) + (... )(a4 - 2as)]s^ 

d, e 

d+f 

Mo oc (...) s + [(...) + (... )(d + |)]s2 

Low-E WlWl s-wave 

Low-E WlWl d-wave 

Coupled-channel A 

e 

M2 oc [(...) + (... )e]s^ 

9 

9 

To, Taoc [{...)+ {...)g]s^ 


Elastic hh-fo 


TABLE II: Relevant combination of the free parameters a, b, a 4 , < 25 , d, e and some useful reactions to extract them from the 
lowest order terms (i.e., s and s^) in a derivative expansion (as well as a few sel ected resonances with the appropriate quantum 
numbers for each channel). The numeric coefficients can be found in appendix A3 so we gloss them over with ellipsis. 


Appendix A: Further details on the effective Lagrangian and 2-body scattering amplitudes 

1. Computation of the amplitudes 

From the Lagrangian in Eq. the following tree-level elastic ujuj ujuj amplitude results, 

= (1 - ^ [2a5S^ + ai{t^ + • 


(Al) 


The one-loop part computation, rather lengthy because of the number of Feynman diagrams, was automated (refs. 
EU), carried out, and reported in m- We obtained 


Aoop(«>^’«) = 3g(4^)2^4 -I- (a^ - l)‘\g{s,t,u)t^ + g{s,u,t)u^)] 

with auxiliary functions 

f{s,t,u) := [20 - 40a^ + 56a‘^ - 72aH + 36b^] 

-h [12 - 24:a^ + SOa'^ - 36aH + 18b^]N, 


1 


-I- [-18 -f 36a^ - 36a^ -f 36a^b - 18b^] log 
-h 3(a^ - if log (+ log(-^ 




gis,t,u) := 26 + 12A£ —91og- 5 - —3 log 

L ^ \ 

where in dimensional regularization D = A — e the pole is contained in 


u 

7^2 


(A 2 ) 


(A3) 

(A4) 


Ag = —h log 47 r — 7 . 
e 

These results coincide with earlier published ones [ 6 ] taking the limit of vanishing light scalar mass there. 
For the ujuj hh amplitude we find, in analogy with Eq. © and at tree level, 

M^vLis, t, u) M 7 e(s> - b)-A + ‘^s^ + ^ (i^ + u^) 


that takes a one-loop correction: 


, .( 1 ) / X o? — h 


yZ yZ 


(A5) 


(A 6 ) 


(A7) 
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where 


f{s,t,u) = -8[-9 + lla^ - 2b] - QN^[-Q + 7a^ - b] 


+ 36(a^ - 1) log-^ 


+ 3(a^ — b) 



+ log 



(AS) 


and the function g is as defined in Eq. (A4). 

Finally, the hh hh elastic amplitude is, at tree-level and keeping only the operator necessary to renormalize the 
one-loop part. 


while the one-loop piece may be written in terms of only one function 

r(s) = 2 + AT, - log 
as 

[ns)s'^ + T{t)e + T{u)u^] . 


(A9) 


(AlO) 


(All) 


2. Renormalization of the amplitudes 


Comparing the tree-level amplitudes in Eqs. ® with the loop ones in Eqs. ( |A2| ), ( |A7| ), ( |A11[ ), we see 

that the divergences in the one-loop pieces can be absorbed just by redefining the couplings a 4 , as, d and e from 
the NLO tree-level Lagrangian. Therefore no renormalizations of a, 6, wave-functions nor (vanishing) masses are 
needed to obtain finite amplitudes (an advantage of dimensional regularization). Our amplitudes are quoted in the 
MS scheme, and the renormalized couplings are 


a. = 


ac — 


9' 


I /I 2 \ 2 

+ 75g^2 ^ ~ ~ 

’+64i5<“ 


'6 + 3&2) 


(A12) 


As a simple limit, the MSM (a = 6 = 1) is renormalizable without any of these additional five couplings (we see that 
they are unnecessary in this case). The case of the Higgsless EWChL corresponds to a = 6 = 0 and then g, d and e 
do not need any renormalization. We also reproduce the well known results for the constants a 4 and m- In more 
generality, the renormalization of a 4 and as agrees with [6]. 

The elastic WBGB amplitude reads, in terms of these renormalized couplings 


A{s,t,u) = ^{1 - a'^) + ^[2al{iJ.)s'^+ +u^)] 

+ t 4 - lOa^ - ISa^b + 96^ + 5)s^ + ^(“^ “ l)^(t^ + u^) 

lOTT^V^ \9 18 

— i (2a^ — 2a^ — 2a^b + b^ + l)s^ log 

+ ^(1-«^)^(s^-u^)log^ 


(A13) 
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While the inelastic ujuj hh amplitude is 


+ 


a^-b . 2d^iii) 2 , e’-(M),,2 


-s + 


-s" + 


+ U^) 


5767r2v4 


—s 


72 - 88a^ + 166 + 36(a^ - 1) log ^ 


+ 3(a^ - 6) log ^ + log ^ 

(a^ _ 5) j 26 - 9 log ^ - 3 log ) 6^ 
(a^ - 6) ( 26 - 9 log - 3 log ^ 


and finally the hh —> hh amplitude may be written as 


T{s,t,u) = 




(A14) 


(A15) 


3(a^ - 6)2 
327r^t’^ 


2 ( 5 ^ + log ^ log ^ log ^ 


Apparently, Eqs. (A13), (A14) and (A15) depend on the renormalization scale /i through the logarithmic terms. But 
they also depend on this arbitrary /i through the renormalized couplings a 4 ... e. 

However, in the absence of wave or mass renormalization, the amplitudes must be observable, and hence /r- 
independent; then we may require that their total derivatives with respect to log vanish. Integrating the resulting 
(very simple) differential equations, we find the renormalization-group evolution equations for the different couplings 
that allow to change the scale 


<(m) 

alih) 

/(m) 

dr{h) 

e’-(M) 


<(Mo) - 

®5(Mo) — 




j^[V-6)^ + 2(l-a“)^]logt 




e(Mo) - 


1 


1927r2 
1 


{o? — b) [{a^ — b) — 6(1 — a^)] log ^ 

/^o 




(A16) 


These equations are diagonal, so the scale evolution does not mix the couplings at NLO in perturbation theory. The 
/i-in varian ce of all the amplitudes has been checked by substituting the /i-evolution of the renormalized couplings in 
Eq. (Alb) into their explicit expressions. 

Erom a practical point of view, we have adopted the values of a^,..., to be as quoted for each example in the 
manuscript at a scale of /i = 3TeV. The dependence on /i is shown in fig. (28) for the / = J = 0 case, and seen to 
be rather moderate. Indeed for a = 0.95, the prefactor of the first Eq. (Alb), say, is 5 x 10“^, so that the scale 
dependence is small. 


3. Detailed partial waves 


E valuat ing the partial-wave projection integral in Eq. ( |13| ) by substituting the renormalized amplitude obtained in 

/?, 


Eq. (A13) for ujuj ujuj provides us with K, D, E constants and H(/i) functions. 

Eor the scalar-isoscalar channel with IJ = 00, the results of m, in terms of the coefficients a, 6, v instead of a 
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jLv=2TeV - jLv=3TeV -------- jLv=4TeV 



FIG. 28: The dependence on the renormalization scale is absorbed throughout in the NLO coefficients. But instead of varying 
them at fixed scale, we can also take the coefficients as fixed (here = 1, 5 = 2 and all the other NLO parameters set to 
zero) and show the dependence on the election of fi. We take for this example the absolute value of the isoscalar amplitude 
(/ = J = 0). There is no qualitative difference in adopting one or another scale. So we have used /x = 3 TeV throughout the 
paper. 


167rv2 


U 


BM = + 768{7a4(Ai) + 


For 


= -460^7(l-«V + 3(«»-i.n 


the vector isovector JJ = 11 amplitude. 


— — 75(a^ — 6)^ + 4608{a4(/i) — 2a5(/x)}7r^ 


1105927r3v4 

' (l-aT- 


For 


E = 

“ 92167r3v4 

the scalar isotensor JJ = 20: 

' 41 - 


■ a 


327rt>2 

3072{2a4(/i) + 


- 28(a^ - b) 


and for the tensor isoscalar IJ = 02, 

Kq2 = 0 

1 


= 1843^,4 + 

‘ (1 - f 


-TT 


__ _ _ 

10247r^x;^ 


Bo2{^i) = 92160^3^4 [320(1 - + 77{a^ - bf + 15360{2a4(Ai) + 

= “460807r3r;4 [10(1-«')"+3(a2 - 6)^] 

Eo2 = 0 . 


(A17) 


(A18) 


(A19) 


(A20) 
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Next we quote a new calculation of the tensor-isotensor / = J = 2 partial wave, that to our knowledge has not 
been reported in the literature. 


K 22 

D 22 

E 22 


9216007r3i;4 

1 


[71(1 - + 77{a^ - bf + 7680{a4(M) + 2a5(Ai)}7r^' 


460807r3v4 


[4(1 -a2)2 +3(a2 -6)2' 


= 0 . 


(A21) 


This exhausts the list of elastic partial waves that are non-vanishing at NLO in perturbation theory, since those with 
angular momentum J = 3 and higher start at 0{s^) and are NNLO in the derivative counting. Needless to say, they 
would be tiny at LHC energies. 

We now give the equivalent results for the inelastic channel-coupling: ujuj hh, with partial waves Mj, starting 
by the scalar one. 


K 

B'M 


D' 


E' 


x/3 

327ru^ 


(a2 


h) 


73 


d(/i) -h 


e(/^) 

3 


A/3(a^ — 6)^ 
92167r^u^ 

V3(a"-6)(l-a") 

5127r^u^ 


. 2 

184327r^u^ 


h) [72(1 - a^) + (a^ 


^)] 


while for the tensor M 2 channel 


K '2 

D'2 

E'2 


0 

e{ii) 83(a^ — 6)^ 

IGOa/Sttu^ ^ 307200V37r3u4 
(g" - hf 
7680A/37r^u^ 

0 . 


(A22) 


(A23) 


At last we quote the elastic hh hh channel amplitude, 
constants 


The To ( 5 ) scalar partial-wave is given by the set of 


while the tensor T 2 requires 


B'iilA 


= 0 


lOfl(M) , (7 - bf 

967rv4 967r3f4 


(a2 - bf 

5127r^v4 

3(7 - bf 
10247r37 


X" = 0 

p/// X ^ g(M) 77(7 - bf 

2 2407ri;4 3072007i;4 

^ (7-6)2 

2 512077 

^2=0. 


(A24) 


(A25) 


By using the evolution equations it is possible to check that all the obtained partial waves are ji independent. 
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Appendix B: Coupled-channel Inverse Amplitude Method 

In this appendix we show how to extend the lAM method to the two-body coupled-channel problem when all 
the partiele species in the various channels are massless. Otherwise this cannot be done because of the presence of 
overlapping cuts. This is well known to happen in the tttt and KK system (see for example [48]) where the KK ^ KK 
left cut terminates at s = 4M|^ — 4M^, which is beyond the tttt ^ tttt threshold branching point located at s = 4M^ 
where the tttt RC starts. Thus, the two cuts overlap. 

In effective theories we can develop the coupled reaction matrix F{s) = according to the chiral/derivative 

expansion 


Fi^{s)=F^\s) + F^f^{s) + ... (Bl) 

where i and j are channel subindices (e.g. j = cjcj, hh ,..) but we have omitted the isospin and momentum indices 
(/, J). As the interactions are assumed to be time-reversal invariant ji. Also in the physical region, i.e. on 

the RC we have: 


ImFij = J2Fu{Fijy 


(B2) 


since coupled-channel unitarity requires the imaginary part of a generic partial wave to receive contributions from all 
allowed intermediate channels. This equation can be written in slightly more compact form as: ImF = FF^ = F^F. 
Now, by using ImF“^ = ImFF“^, the unitarity condition on the RC can be written as: 

ImF-^ = -l. (B3) 

However in the effective theory this condition is only satisfied perturbatively, so that at one-loop precision 

ImF^P {F^py . (B4) 


As we saw in Eq. (15), the lowest-/J NLO partial-waves take the general form 


Fh 


NLO 


(s) = F^^is) + F^pis) = Ki^s + + Ai log ■ 


SO that the perturbative unitarity of Eq. (B4) on the physical RC requires 

Eij = . 


Eij log 


(B5) 


(B6) 


Now, by following the same steps as in the single channel case in section we can obtain a twice-subtracted DR for 
the Ffj^^{s). Next we introduce the inverse amplitude matrix function as: 


w = 


(B7) 


The essential point here is that, as all the particles are massless, the analytical structure of all the matrix elements 
Fij{s) is the same, namely a LC and a RC starting at the origin. This structure is also shared by each of the F^^{s) 
and Wij{s) matrix elements. Had the masses of the particles appearing in the various channels, and consequently the 
cut structure, been different, the Wij{s) matrix elements would mix and possibly overlap the different left and right 
cuts. This would produce spurious contributions to the imaginary part of the partial-waves in the physical region. 
Considering only massles s par ticles ensures that we will not have this kind of spurious contributions. Extracting the 
imaginary part from Eq. (B7) on the RC we obtain: 


ImW = Im = -F‘-°'>F‘~°^ = - Im 


(B8) 


where we have used Eq. (B4). Then we have: 


ImWijis) = -lmFlp(s) = -E-. 


(B9) 
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By using it in a twice-subtracted DR for each Wij{s) matrix element, and assuming that no poles appear when 
inverting the F{s) matrix, we obtain on the LC at NLO: 

ImVF --lmF(i), (BIO) 

as we did in the single channel case. Thus we finally find: 


—s 


W,,{s) ^ + A, log ^ + A, log — ) = - F^;\s) 






p(i)/ 


■ 


and then we arrive to the lAM formula for massless coupled channels in Eq. (45): 

-1 


^lAM ^ 


^(0) F^^\ 


(Bll) 


(B12) 


As discussed already in subsection IIIB this matrix is exactly unitary, i.e. ImF^^^ = ^iam ^^lAM^t 
it is also compatible with the NLO approximation as 


Fl^^{s)=Fl°\s) + F^^p{s) + 0{s% 


(B13) 


In addition all the elements have the same proper analytical structure (left and right cuts). This makes the analytical 
continuation to the second Riemann sheet of the different amplitudes possible, and eventually the presence of poles 
there, that could be understood as dynamical resonances for some regions of the parameter space. 

Once again, this construction is possible only because all the particles are assumed to be massless. This is a good 
approximation because we are using the ET and the Landau gauge in our computations and also because the physical 
h mass — 125 GeV is close to Mw and M^. However, if we had taken different from zero the hh hh channel 
would have had a LC ending at the positive value 4M|. In this situation the lAM multichannel method produces a 
spurious imaginary part in the physical region of the channel ujuj ujuj ranging from s = 0 to s = 4M| thus spoiling 
unitarity in that region. 


Appendix C: The N/D solution 


In this appendix we detail some computations needed to construct an approximation to the N/D method in sec¬ 
tion IV The dispersion relation for the denominator A in Eq. (48), under the approximation N{){s) + 


results in 


Dq[s) — 1 + his + h2S^ -( KI 2 

TT 


D 




D^E 

where the IR and UV regularized In and A integrals are respectively defined as 


and 


In{s;m,A) = / 

J m 

I'n{s]m,A,ii) = / 
J m 


ds' 


2 — s — ie) 


^ ds' log ^ 

2 — s — ie) 


(Cl) 


(C2) 


(C3) 


These integrals can be computed and are UV convergent thanks to n being positive because of the subtractions in 
the dispersion relation. Thus, taking A ^ 00 directly, one finds 


12 ( 5 ; m,oc) = 

1 

r s / 

9 

— + log 1 


Im^ V 

/i(s;m, oc) = 

1 

/ \ 

— log 1 -^ 

s 

\ / 

A(«s; m, 00 ,/i) = 

1 

s 





m" 




-log ^ ^ - log ^ log 1- j - Li 2 ( — + A log ^ 


m" 


m" 


(C4) 
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where 1/^2 (^) is the dilogarithm function. Therefore we can write: 


yl(s) - An/d(s) ^ +^l(s) - 

^ l + h^s + h2s‘^-^T{s) 


(C5) 


where: 


T{s) = Kh + p^Dh + DI[. 


Now it is not difficult to check that for small m 


Do{s) = l + his+ ——— log f + O(s^) 
TT V } 


(C6) 


(C7) 


SO that matching the dispersion relation to perturbation theory sets the hi subtraction constant, the correct choice 
being 


hi = hi{m) = 


K 


7r{D + E) 


B{m) 


so we have 




(C8) 


(C9) 


and then 


^n/d(<^) 


No{s) 

Do{s) 


+ ^l{s) + Aji{s) + O(s^) 


(CIO) 


which reproduces the NLO computation, 
show that for small enough IR cutoff m: 


From the integrals above in Eq. (C2) and (C3) it is also not difficult to 


ImT(5) 




(Cll) 


on the RC so that the ^n/d(<^) partial waves fulffil exact elastic unitarity. 

In order to have a clearer mathematical description of the amplitude obtained it is useful to introduce an additional 
subtraction constant. Thus we define 


II = H{m) = h2{m)7T +E+ . 


(C12) 


Then it is very easy to show that 


^ 0 ( 5 ) = 1 


Ar{s) 

A(^){s) 


H{m) 


D 


D 


D^E 




Bill) log ^ — log ^ - log ^ 


/i" 


11^ 


(C13) 


As usual H = Him) must be considered a renormalized parameter at the scale m. By demanding D^is) to be 
independent of this scale we find the renormalization equation: 


m 


.dHim) 


D 


-Bim) 


(C14) 


dm‘^ D ^ E 

which upon integration leads to an evolution equation characteristic of an NNLO parameter in perturbation theory, 


(m) = /7(mo) + log ^ log2 ^ . 

U ^ IIq Z IIq 

Finally we can remove the IR cutoff m from D^is) to find the /i independent equation: 


o 2 r 


^ + 7 


D , —s D ^ 2 —s 

— 5(M)log^ + -log 


(C15) 


(C16) 
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The first term, 1 here, corresponds to LO in perturbation theory (the s power being contained in the numerator N). 
The Ar term contains the NLO physics, and finally the method has generated an NNLO piece that is necessary to 
have the correct analytic properties. Thus the renormalized constant can contain contributions from the NNLO 

chiral couplings. However it is possible to neglect these contributions in a consistent way by choosing: 



(C17) 


which, at it is easy to check, satisfies the above evolution equation. With this choice the partial wave denominator 
takes the simpler form: 




(C18) 


which we have used in the main text [Eq. ([60|)]. By using the the Al{s) and Ar{s) definitions this denominator can 
also be written as: 


Do{s) = 1 - 


Ar{s) Al{-s)Ar{s) 


H(o)(5) 2 ( H ( o ))2 

Notice that this denominator has not any LC as it must be the case. 


(C19) 


Appendix D: Numerical extraction of coupled-channel poles in complex s 


Here we describe very briefly the numeric finding of resonances and their parameter extraction as poles of the 
amplitude in the complex-s plane, and also to assess when violations of causality occur (finding poles in the first 
Riemann sheet instead of the second). 

An accurate method is the use of a Cauchy line integral around a closed path; a finite value indicates that some 
pole has been enclosed. The difficulty comes from having two coupled-channels, though it is not severe since we have 
taken all particles as massless so that the cuts of the two channels start at the same point (the origin in the complex 
s plane). 

We find convenient to use an integration contour shaped as a half-circle out at a radius R = (3TeV)^ (roughly the 
range of validity of the unitarization methods considered in this work), closed by a segment of the imaginary axis, 
and parametrized in terms of a dummy integration variable t. 


tW = 



for t G [0,1] 
for t G (1, 2] 


(Dl) 


Cauchy’s theorem states that if a function A{s) has N poles on points Si (i = 1,..., A) within the region enclosed by 
7 , the value of the line integrals 



N 

i=l 


(D2) 


is given by the respective pole residues Ho(si) (there is no reason to expect double poles in our NLO-based compu¬ 
tation). The A{s) function generically stands for any of the considered scattering amplitudes. Since the low-energy 
perturbative interactions are weak, we do not look for bound states and thus restrict Res > 0. For 1ms > 0, the 
logarithms in A{s) are evaluated on the first Riemann sheet; for 1ms < 0, on the second. 

Now, the circuit in Eq. |D2| is taken on the second Riemann sheet. Thus, it immediately captures all poles on its 
lower half-plane and also on the upper half-plane (common to both first and second sheets). Poles in the lower half 
plane of the first Riemann sheet are outside the contour. Still, we can detect them because they occur simultaneously 
with a pole on the upper half plane, as we now argue. 

Because of the analytical properties of scattering amplitudes Aj{s) on the first Riemann sheet (analyticity on the 
upper half plane plus cut along the positive real-s axis), Schwarz reflection applies and [A/( 5 *)]* = Ai{s). So, every 
pole on the first Riemann sheet below the real axis (Im s < 0) implies the presence of a pole at s' = s* over the real 
axis (Im s' > 0). 

In contrast, by definition, in the second Riemann sheet A{s) will be analytic on 1ms = 0, Res > 0, so the pole on 
the lower-half plane of the second sheet does not reflect on the upper half-plane. 
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Thus, the path given by Eq. is sufficient to detect all poles generated with the lAM method within its range 
of validity, on both the first and second Riemann sheets (those with Im s < 0 are on the second, those with Im 5 > 0 
tag a pair on the first sheet, respectively). 

For each studied parameter set, three integrals (Ik) of the family in Eq. 


D2 


have been computed, with /c = 0, 1, 2. 
If no pole lies inside the contour, the value of all these integrals is zero. Next, if we have only one pole at position s, 
we can equate two ratios of these Ik integrals. 


. h .2 h 

s=Y.s =- 


/i = 10 I 2 


(D3) 


For a larger contained-pole count, N > 1 {e.g. one on the first and one on the second Riemann sheets), it wouldn’t be 
generally true that = / 0 / 2 , so we can use this relation as a check of whether there is exactly one pole there. Should 
it fail, a more detailed study would be necessary. However, it would still be possible to compute the position of an 
arbitrary (but finite) number of poles by computing integrals of increasing order and solving the non-linear equation 
system Ik = 27ri ^ s^Ao{si). For the particular case of = 2, the expressions would still be analytic. And, in 
particular. 


51,2 


ihh - loh) ± VOih - hh? - 4(/? - hh)m - hh) 
2{ll - hh) 


We also record here some analytical expressions to find the location of 3 poles from the h, k = 0,. 
It is best quoted in terms of several auxiliary quantities, namely 


A 

A-A 
A-B 
A-C 
A 




2hhh - hit - Ilh + hhh 


= -Ilh + hh + hhh - hhh - hh + / 1 / 3/6 
= -hh + hh+hhh - hh - hhh+hhh 


— J-A 


2hhh - hh - hii 


= -A^B^ +AB-^ + 4:A^C 


r = 


-2A-^ + 9AB - 27C + 3V3A 


i-'s 

3-^6 + hhh 
ISisC + 27(7^ 
1/3 


,., 5 integrals. 

(D4) 

(D5) 

(D6) 

(D7) 

(D8) 

(D9) 


Then the pole locations become 


A 21/3(3.5 - i^) r 

"3 3f 3 ■ 21/3 

A 3{1 + iy^){B + A^) (l-iV3)r 

"3 3 ■ 22/3r 6 ■ 21/3 

A 3{1 - iV^){B + A^) (l + i\/3)r 
s’ 3 • 22/3r 6 • 21/3 


(DIO) 

(Dll) 
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